IGPG-06/03-2 

Quantum Nature of the Big Bang: 
An Analytical and Numerical Investigation 

Abhay Ashtekar^'^'^jj Tomasz Pawlowski^ tJ and Parampreet Singh^'xl 

^Institute for Gravitational Physics and Geometry 

Physics Department, Penn State, University Park, PA 16802, U.S.A. 

"^Inter-University Gentre for Astronomy and Astrophysics 

post hag 4, Ganeshkhind, Pune 4II 017, India 

^ Isaac Newton Institute for Mathematical Sciences, 

^ : 20 Glarkson Road, Gamhridge GB3 OEH, UK 

O ■ Abstract 

Analytical and numerical methods are developed to analyze the quantum nature of the big bang 

^ . in the setting of loop quantum cosmology. They enable one to explore the effects of quantum 

geometry both on the gravitational and matter sectors and significantly extend the known results 

On I on the resolution of the big bang singularity. Specifically, the following results are established for 

the homogeneous isotropic model with a massless scalar field: i) the scalar field is shown to serve as 

cn ' an internal clock, thereby providing a detailed realization of the 'emergent time' idea; ii) the phys- 

jK^ . ical Hilbert space, Dirac observables and semi-classical states are constructed rigorously; iii) the 
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Hamiltonian constraint is solved numerically to show that the big bang is replaced by a big bounce. 
Thanks to the non-perturbative, background independent methods, unlike in other approaches the 
r^ ' quantum evolution is deterministic across the deep Planck regime. Our constructions also provide 
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a conceptual framework and technical tools which can be used in more general models. In this 
sense, they provide foundations for analyzing physical issues associated with the Planck regime of 



qh, loop quantum cosmology as a whole 
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I. INTRODUCTION 

Loop quantum gravity (LQG) is a background independent, non-perturbative approach 
to quantum gravity [l|, 12; |3|. It is therefore well-suited for the analysis of certain long 
standing questions on the quantum nature of the big bang. Examples of such questions are: 

• How close to the Big Bang does a smooth space-time of general relativity make sense? In 
particular, can one show from first principles that this approximation is valid at the onset 
of inflation? 

• Is the Big-Bang singularity naturally resolved by quantum gravity? Or, is some external 
input such as a new principle or a boundary condition at the Big Bang essential? 

• Is the quantum evolution across the 'singularity' deterministic? Since one needs a fully 
non-perturbative framework to answer this question in the affirmative, in the Pre-Big-Bang 
1^ and Ekpyrotic /Cyclic j5l, |6i] scenarios, for example, so far the answer has been in the 
negative |7|, |lO| . 

• If the singularity is resolved, what is on the 'other side'? Is there just a 'quantum foam', 
far removed from any classical space-time (as suggested, e.g., in [ll|), or, is there another 
large, classical universe (as suggested, e.g., in jj,|5|,|6|)? 

Over the years, these and related issues had been generally relegated to the 'wish list' of 
what one would like the future, satisfactory quantum gravity theory to eventually address. 
However, over the past five years, thanks to the seminal ideas introduced by Bojowald and 
others, notable progress was made on such questions in the context of symmetry reduced, 
minisuperspaces |l2| • In particular, it was found that Riemannian quantum geometry, which 
comes on its own in the Planck regime, has just the right features to resolve the big-bang sin- 
gularity jl3| in a precise manner 1J| . However, the physical ramifications of this resolution 



— in particular the answer to what is 'on the other side' — have not been worked out. It is 
therefore natural to ask whether one can complete that analysis and systematically address 
the questions listed above, at least in the limited context of simple cosmological models. It 
turns out that the answer is in fact in the affirmative. A brief summary of arguments leading 
to this conclusion appeared in [15[. The purpose of this paper is to provide the detailed 
constructions and numerical simulations that underlie those results. 

Let us begin with a brief summary of the main results of loop quantum cosmology (LQG). 
(For a comprehensive survey, see, e.g., |l3|-) They can be divided in two broad classes: 
i) singularity resolution based on exact quantum equations (see e.g. |l3|, [ij, |l6|, |l7|, |l8|, |l9|, 
HEI), and. 



ii) phenomenological predictions based on 'effective' equations (see e.g. [2j, |23, U3, |25|, |2( 

3,111 m, H H H H in ) . 



Results in the first category make a crucial use of the effects of quantum geometry on 
the gravitational part of the Hamiltonian constraint. Because of these effects, the quantum 
'evolution' is now dictated by a second order difference equation rather than the second 
order differential equation of the Wheeler-DeWitt (WDW) theory. Nonetheless, the intuitive 
idea of regarding the scale factor a as 'internal time' was maintained. The difference is 
that the 'evolution' now occurs in discrete steps. As explained in detail in Sec. IIIB[ this 
discreteness descends directly from the quantum nature of geometry in LQG; in particular, 
the step size is dictated by the lowest non-zero eigenvalue of the area operator. When the 
universe is large, the WDW differential equation is an excellent approximation to the 'more 
fundamental' difference equation. However, in the Planck regime, there are major deviations. 
In particular, while the classical singularity generically persists in the WDW theory without 



additional inputs, this is not the case in LQG. This difference does not arise simply because 
the discrete 'evolution' enables one to 'jump over' the classical singularity. Indeed, even 
when the discrete evolution passes through the point a = 0, the difference equation remains 
well-defined and enables one to 'evolve' any initial data across this classically singular point. 
It is in this sense that the singularity was said to be resolved. 

In the second category of results, by and large the focus was on quantum geometry 
modifications of the matter part of the Hamiltonian constraint. The main idea is to work 
with approximation schemes which encode the idea of semi-classicality and to incorporate 
quantum geometry effects on the matter Hamiltonian by adding suitable 'effective' terms to 
the classical Hamiltonian constraint. The hope is that dynamics generated by the 'quantum- 
corrected' classical Hamiltonian would have a significant domain of validity and provide a 
physical understanding of certain aspects of the full quantum evolution. This strategy has 
led to a number of interesting insights. For example, using effective equations, it was argued 
that the singularity resolution is generic for all homogeneous models; that there is an early 
phase of inflation driven by quantum gravity effects; that this phase leads to a reduction of 
power spectrum in CMB at large angular scales; and that quantum geometry effects suppress 
the classical chaotic behavior of Bianchi IX models near classical singularity. 

Attractive as these results are, important limitations have persisted. Let us begin with 
the results on singularity resolution. As in most non-trivially constrained systems, solutions 
to the quantum constraint fail to be normalizable in the kinematical Hilbert space on which 
the constraint operators are well-defined. Therefore one has to endow physical states (i.e. 
solutions to the Hamiltonian constraint) with a new, physical inner product. Systematic 
strategies — e.g., the powerful 'group averaging procedure' J35|, |36| — typically require that 



the constraint be represented by a self-adjoint operator on the kinematic Hilbert space while, 
so far, most of the detailed discussions of singularity resolution are based on Hamiltonian 
constraints which do not have this property. Consequently, the space of solutions was not 
endowed with a Hilbert space structure. This in turn meant that one could not introduce 
Dirac observables nor physical semi-classical states. Therefore, as pointed out, e.g. in J37[, 
the physical meaning of singularity resolution remained somewhat obscure. In particular, 
even in simple models, there was no clear-cut answer as to what the universe did 'before' the 
big-bang. Was there a genuine 'quantum foam' or was the quantum state peaked at a large 
classical universe on the 'other side' ? In absence of a physical Hilbert space, the second and 
the third questions posed in the beginning of this section could only be answered partially 
and the first and the fourth remained unanswered. 

The phenomenological predictions have physically interesting features and serve as valu- 
able guides for future research. However, the current form of this analysis also faces some 
important limitations. Many of these discussions focus only on the quantum geometry mod- 
ifications of the matter Hamiltonian. This strategy has provided new insights and does serve 
as a useful starting point. However, there is no a priori reason to believe that it is consis- 
tent to ignore modifications of the gravitational part of the Hamiltonian and retain only 
the matter modifications. Conclusions drawn from such analysis can be taken as attractive 
suggestions, calling for more careful investigations, rather than firm predictions of LQC, to 
be compared with observations. On the conceptual side, a number of semi-classical approx- 
imations are made while deriving the effective equations. Many of them are largely violated 
in the Planck regime. Therefore it is difficult to regard conclusions drawn from effective 
equations on singularity avoidance, and on the fate of the universe beyond, as reliable. It 
does happen surprisingly often in physics that approximation schemes turn out to work even 



outside the regimes for which they were originally intended. But there is no a priori reason 
to think that this must happen. To develop intuition on the validity of approximations, it 
is essential to make, in at least a few simple models, detailed comparisons of predictions of 
effective equations with those of quantum equations that are being approximated. To sum- 
marize, while LQC has led to significant progress by opening new avenues and by indicating 
how qualitatively new and physically desirable results can arise, the program has remained 
incomplete even within the realm of symmetry reduced, mini-superspace models. 

The goal of this paper is to complete the program in the simplest of models by using a 
combination of analytical and numerical methods. The resulting theory will enable us to 
answer, in the context of these models, the questions raised in the beginning of this section. 
Specifically, we will show from first principles that: i) a classical space-time continuum is 
an excellent approximation till very early times; ii) the singularity is resolved in the sense 
that a complete set of Dirac observables on the physical Hilbert space remains well-defined 
throughout the evolution; iii) the big-bang is replaced by a big-bounce in the quantum theory; 
iv) there is a large classical universe on the 'other side', and, v) the evolution bridging 
the two classical branches is deterministic, thanks to the background independence and 
non-perturbative methods. While the paper is primarily concerned with basic conceptual 
and computational issues, our constructions also provide some tools for a more systematic 
analysis of phenomenological questions. Finally, our approach can be used in more general 
models. In particular, our constructions can be used also for anisotropic models and for 
models in which the scalar field has a potential, although certain conceptual subtleties have 
to be handled carefully and, more importantly, the subsequent numerical analysis is likely 
to be significantly more complicated. Nonetheless, in a rather well defined sense, these 
constructions provide a foundation from which one can systematically analyze the Planck 
regime in LQC well beyond the specific model discussed in detail. 

The main ideas of our analysis can be summarized as follows. First, our Hamiltonian 
constraint is self-adjoint on the kinematical (more precisely, auxiliary) Hilbert space. Second, 
we use the scalar field as 'internal time'. In the classical theory of /c = models with a 
massless scalar field 0, the scale factor a as well as are monotonic functions of time in 
any given solution. In the k = l case, on the other hand, since the universe recollapses, only 
has this property. Therefore, it seems more natural to use as 'internal time', which 
does not refer to space-time coordinates or any other auxiliary structure. It turns out that 
is well-suited to be 'emergent time' also in the quantum theory. Indeed, our self-adjoint 
Hamiltonian constraint is of the form 

where G does not involve and is a positive, self-adjoint, difference operator on the auxiliary 
Hilbert space of quantum geometry. Hence, the quantum Hamiltonian constraint can be 
readily regarded as a means to 'evolve' the wave function with respect to 0. Moreover, this 
interpretation makes the group averaging procedure similar to that used in the quantization 
of a 'free' particle in a static space-time, and therefore conceptually more transparent. Third, 
we can carry out the group averaging and arrive at the physical inner product. Fourth, 
we identify complete sets of Dirac observables on the physical Hilbert space. One such 
observable is provided by the momentum p^ conjugate to the scalar field and a set of 



them by aj,^^, the scale factor at the 'instants of emergent time' 0o- ^ Fifth, we construct 
semi-classical states which are peaked at values of these observables at late times, when 
the universe is large. Finally we 'evolve' them backwards using the Hamiltonian constraint 
using the adaptive step, 4th order Runge-Kutta method. The numerical tools are adequate 
to keep track not only of how the peak of the wave function 'evolves' but also of fluctuations 
in the Dirac observables in the course of 'evolution'. A variety of numerical simulations have 
been performed and the existence of the bounce is robust. 

The paper is organized as follows. In section |n] we summarize the framework of LQC 
in the homogeneous, isotropic setting, keeping the matter field generic. We first summarize 
the kinematical structure 1J| , highlighting the origin of the qualitative difference between 



the WDW theory and LQC already at this level. We then provide a self-contained and 
systematic derivation of the quantum Hamiltonian constraint, spelling out the underlying 
assumptions. In section ITTTl we restrict the matter field to be a massless scalar field, present 
the detailed WDW theory and show that the singularity is not resolved. The choice of a 
massless scalar field in the detailed analysis was motivated by two considerations. First, it 
makes the basic constructions easier to present and the numerical simulations are substan- 
tially simpler. More importantly, in the massless case every classical solution is singular, 
whence the singularity resolution by LQC is perhaps the clearest illustration of the effects of 
quantum geometry. In section HVl we discuss this model in detail within the LQC framework, 
emphasizing the dynamical differences from the WDW theory. Specifically, we present the 
general solution to the quantum constraint, construct the physical Hilbert space and in- 
troduce operators corresponding to Dirac observables. In section |V| we present the results 
of our numerical simulations in some detail. Solutions to the Hamiltonian constraint are 
constructed using two different methods, one using a Fast Fourier Transform and another by 
'evolving' initial data using the adaptive step, 4th order Runge-Kutta method. To further 
ensure the robustness of conclusions, the initial data (at late 'times' when the universe is 
large) is specified in three different ways, reflecting three natural choices in the construction 
of semi-classical states. In all cases, the classical big-bang is replaced by a quantum big- 
bounce and the two 'classical branches' are joined by a deterministic quantum evolution. 
Section IVTl compares and contrasts the main results with those in the literature. 

Issues which are closely related to (but are not an integral part of) the main results 
are discussed in three appendices. Appendix [X] is devoted to certain heuristics on effective 
equations and uncertainty relations which provide a physical intuition for 'mechanisms' 
underlying certain constructions and results. In Appendix |B1 we discuss technical aspects of 
numerical simulations which are important but whose inclusion in the main text would have 
broken the flow of the argument. Finally in Appendix we present an alternate physical 
Hilbert space which can be constructed by exploiting certain special features of LQC which 
are not found in the general setting of constrained systems. This space is more closely related 
to the WDW theory and could exist also in more general contexts. Its existence illustrates 
an alternate way to address semi-classical issues which may well be useful in full LQC. 



^ As we will see in section ^ to construct a Hamiltonian framework in the open model, one has to fix a 
fiducial cell. The scale factor a (and the momentum p conjugate to the gravitational connection introduced 
later) refers to the volume of this cell. Alternatively, one can avoid the reference to the fiducial cell by 
fixing a (f>o and considering the ratios aj^/a^^ as Dirac observables. However, for simplicity of presentation 
we will not follow this route. 



II. LQC IN THE k = HOMOGENEOUS, ISOTROPIC SETTING 

This section is divided into three parts. To make the paper self-contained, in the first we 
provide a brief summary of the kinematical framework, emphasizing the conceptual structure 
that distinguishes LQC from the WDW theory. In the second, we introduce the self-adjoint 
Hamiltonian constraints and their WDW limits. By spelling out the underlying assump- 
tions clearly, we also pave the way for the construction of a more satisfactory Hamiltonian 
constraint in [38|-^ In the third part, we first list issues that must be addressed to extract 
physics from this general program and then spell out the model considered in the rest of the 
paper. 

A. Kinematics 

1. Classical phase space 

In the k = case, the spatial 3-manifold M is topologically M^, endowed with the action 
of the Euclidean group. This action can be used to introduce on M a fiducial, flat metric °qab 
and an associated constant orthonormal triad °e" and co-triad lu*. In full general relativity, 
the gravitational phase space consists of pairs (A^, E!^) of fields on M, where A^ is an SU(2) 
connection and E"' an orthonormal triad (or moving frame) with density weight 1 [1]. As 
one would expect, in the present homogeneous, isotropic context, one can fix the gauge and 
spatial diffeomorphism freedom in such a way that A^ has only one independent component, 
c, and El, only one independent component, p: 

A = c%'Ti, E = p^%iT\ (2.1) 

where c and p are constants and the density weight of E has been absorbed in the determinant 
of the fiducial metric.^ Denote by T^^^^ the subspace of the gravitational phase space Fgrav 
defined by ()2.H) . Because M is non-compact and our fields are spatially homogeneous, 
various integrals featuring in the Hamiltonian framework of the full theory diverge. However, 
the presence of spatial homogeneity enables us to bypass this problem in a natural fashion: 
Fix a 'cell' V adapted to the fiducial triad and restrict all integrations to this cell. (For 
simplicity, we will assume that this cell is cubical with respect to °qah-) Then the gravitational 
symplectic structure figrav on Fgrav and fundamental Poisson brackets are given by 14| : 



f^frav = ^ dg A dp, and {c, P} = ^ (2.2) 



^ Our conventions are somewhat different from those in the Hterature, especially [lj|. First, we follow 
the standard quantum gravity convention and set ^pj = Gh (rather than SirGh). Second, we follow the 
general convention in geometry and set the volume element e on M to be e :— y^\ det-E| (rather than 
e := y^ I det E \ sgn(det E) ) . This gives rise to some differences in factors of sgnp in various terms in the 
expression of the Hamiltonian constraint. Finally, the role of the the minimum non-zero eigenvalue of 
area is spelled out in detail, and the typographical error in the expression of iIq that features in the 
Hamiltonian constraint is corrected. 

^ Our conventions are such that Ti tj — 2^ijkT^ — i^ij- Thus, 2iTk — ak, where ai are the Pauli matrices. 



where Vo is the volume of V with respect to the auxihary metric °qab- Finally, there is a 
freedom to rescale the fiducial metric °qab by a constant and the canonical variables c, p fail 
to be invariant under this rescaling. But one can exploit the availability of the elementary 
cell V to eliminate this additional 'gauge' freedom. For, 

1 2 

c = Vo^c and p = Vo^p (2.3) 

are independent of the choice of the fiducial metric °qab- Using (c, p), the symplectic structure 
and the fundamental Poisson bracket can be expressed as 

^frav = ^^ dc A dp and {c, p} = ^^ . (2.4) 

Since these expressions are now independent of the volume Vo of the cell V and make no 
reference to the fiducial metric °qab, it is natural to regard the pair (c, p) as the basic canonical 
variables on T^^^^. In terms of p, the physical triad and co-triad are given by: 

el = (sgnp)|p|-5 (v) X), and cul = {sgnp)\p\i {Vo''' "o^:) (2.5) 



The function sgnp arises because in connection-dynamics the phase space contains triads 
with both orientations, and since we have fixed a fiducial triad °e", the orientation of the 
physical triad e^ is captured in the sign of p. (As in the full theory, we also allow degenerate 
co-triads which now correspond to p = 0, for which the triad vanishes.) 

Finally, note that although we have introduced an elementary cell V and restricted all 
integrals to this cell, the spatial topology is still M^ and not T^. Had the topology been 
toroidal, connections with non-trivial holonomy around the three circles would have en- 
larged the configuration space and the phase space would then have inherited additional 
components. 

2. Quantum kinematics 

To construct quantum kinematics, one has to select a set of elementary observables which 
are to have unambiguous operator analogs. In non-relativistic quantum mechanics they are 
taken to be x,p. One might first imagine using c, p in their place. This would be analogous 
to the procedure adopted in the WDW theory. However, unlike in geometrodynamics, LQG 
now provides a well-defined kinematical framework for full general relativity, without any 
symmetry reduction. Therefore, in the passage to the quantum theory, we do not wish to 
treat the reduced theory as a system in its own right but follow the procedure used in the 
full theory. There, the elementary variables are: i) holonomies he defined by the connection 
A^ along edges e, and, ii) the fluxes of triads Ef across 2-surfaces S. In the present case, this 
naturally suggests that we use: i) the holonomies h/^ along straight edges (/x °e^) defined by 
the connection A = [c/ \^) to^Ti , and ii) the momentum p itself ^]. Now, the holonomy 
along the kth edge is given by: 

h[^) = cos ^1 + 2 sin ^Tfc (2.6) 

where I is the identity 2X2 matrix. Therefore, the elementary configuration variables can 
be taken to be exp(i/^c/2) of c. These are called almost periodic functions of c because /x is 



an arbitrary real number (positive if the edge is oriented along the fiducial triad vector °ef, 
and negative if it is oriented in the opposite direction). The theory of these functions was 
first analyzed by the mathematician Harold Bohr (who was Niels' brother). 

Thus, in LQC one takes e*^*^/^ and p as the elementary classical variables which are to 
have unambiguous operator analogs. They are closed under the Poisson bracket on T^^^^. 
Therefore, as in quantum mechanics, one can construct an abstract ^ir-algebra a generated 
by e*'^^/^ and p, subject to the canonical commutation relations. The main task in quantum 
kinematics is to find the appropriate representation of a. 

In this task, one again follows the full theory. There, surprisingly powerful theorems are 
now available: By appealing to the background independence of the theory, one can select 
an irreducible representation of the holonomy-fiux algebra uniquely J39|, |40|. The unique 



representation was constructed already in the mid-nineties and is therefore well understood 
36l l4lL l42l . I43I I44I l4q . In this representation, there are well defined operators corresponding 



to the triad fiuxes and holonomies, but the connection itself does not lead to a well-defined 
operator. Since one follows the full theory in LQC, the resulting representation of a also 
has well-defined operators corresponding to p and almost periodic functions of c (and hence, 
holonomies), but there is no operator corresponding to c itself. 

The Hilbert space underlying this representation is H^^'^ = L^(MBohr, c^AtBohr), where 
l^Bohr is the Bohr compactification of the real line and /iBohr is the Haar measure on it. It 
can be characterized as the Cauchy completion of the space of continuous functions /(c) on 
the real line with finite norm, defined by: 

WfW = lim ^ r ffdc. (2.7) 

An orthonormal basis in Ti^in^ is given by the almost periodic functions of the connection, 
-^(m) ('^X-i," e^^'^l'^ . (The N(^^){c) are the LQC analogs of the spin network functions in full 
LQG BliJ). They satisfy the relation 

(iV(^)liV(^o) = {e-^\e-^) = 5,,,'. (2.8) 

Note that, although the basis is of the plane wave type, the right side has a Kronecker delta, 
rather than the Dirac distribution. Therefore a generic element \l/ of 'H^^ can be expanded 
as a countable sum \l/(c) = Ylik ^kN^^,^) where the complex coefficients a^ are subject to 
Sfc I'^fcP < ^^- Consequently, the intersection between TC^^^ and the more familiar Hilbert 
space L^(R, (ic) of quantum mechanics (or of the WDW theory) consists only of the zero 
function. Thus, already at the kinematic level, the LQC Hilbert space is very different from 
that used in the WDW theory. 

The action of the fundamental operators, however, is the familiar one. The configuration 
operator acts by multiplication and the momentum by differentiation: 

iV(^)^(c) = exp^^(c), and p^(c) = -z^-^ — (2.9) 

where, as usual, £pi = Gh. The first of these provides a 1-parameter family of unitary 
operators on TC^^^ while the second is self-adjoint. 

It is often convenient to use the Dirac bra-ket notation and set e'^^^^'^ = (c|yu). In this 
notation, the eigenstates of p are simply the basis vectors |/x): 

p|/.) = ^/i|/i). (2.10) 
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Finally, since the operator V representing the volume of the elementary cell V is given by 
V = |]3|2, the basis vectors are also eigenstates of V: 



\/|/.) = (^-^l/^IJ 4il/i)- (2.11) 

The algebra a, of course, also admits the familiar representation on L^(M, dc). Indeed, 
as we will see in section UTTl this is precisely the 'Schrodinger representation' underlying the 
WDW theory. The LQC representation outlined above is unitarily inequivalent. This may 
seem surprising at first in the light of the von-Neumann uniqueness theorem of quantum 
mechanics. The LQG representation evades that theorem because there is no operator c 
corresponding to the connection component c itself. Put differently, the theorem requires 
that the unitary operators iV(^) be weakly continuous in fi, while our operators on H^j^ are 
not. (For further discussion, see |48[.) 



B. The Hamiltonian constraint 

1 . Strategy 

Because of spatial flatness, the gravitational part of the Hamiltonian constraint of full 
general relativity simplifies and assumes the form: 

Cgrav = -7"' / d'xNe,,kFl, e-'E^'^E'' (2.12) 

Jv 

where e := a/| det E\, and where we have restricted the integral to our elementary cell V. 
Because of spatial homogeneity the lapse N is constant and we will set it to one. 

To obtain the corresponding constraint operator, we need to first express the integrand 
in terms of our elementary phase space functions he,p and their Poisson brackets. The term 
involving triad can be treated using the Thiemann strategy [3|, |5l| : 

8,,, e-' E^^E'' = V ^'^^^^ , %»^^ "u', Tr (/^^^ {/.^^-\ V} r,) (2.13) 

where tvl! is the holonomy along the edge parallel to the fcth basis vector of length [i^yo 
with respect to °qah-, and V = \pf^'^ is the volume function on the phase space. While the 
right side of ()2.13p involves /Xq, it provides an exact expression for the left side which is 
independent of the value of /Iq. 

For the field strength F^^, we use the standard strategy used in gauge theories. Consider 
a square Djj in the i-j plane spanned by two of the triad vectors °e^, each of whose sides 

1/3 

has length HoVo with respect to the fiducial metric "qat- Then, 'the ab component' of the 
curvature is given by: 

F^, = -2 hm^ Tr -^^ ] r^%\ ^^i (2.14) 



Mo- , „„„ 



t^lVo' 



where the holonomy h^° around the square Djj is just the product of holonomies ()2.6|) 
along the four edges of Dj^: 

/i2^°) = h[''°^hf°\h[''"^)-\h^^"^)-^ . (2.15) 

By adding the two terms and simphfying, Cg^av can be expressed as: 

Cgrav = hm C^J^^l where 



Mo 



4(sgnp) Y^i,-fcrp^. ,(/..), (Mo), (/^o)-1,(a'c)-1, (Mo) r,(Mo)-l 



ci^^l = -PT57^E^^''Tr(/^r^C^^r^"'^r^"'^r^{/^r^"\n) (2.16) 

Since Cgmv is now expressed entirely in terms of holonomies and p and their Poisson bracket, 
it is straightforward to write the corresponding quantum operator on H^^^: 

24z(sgnp) .2 /. f^oCy floC fioCy . fioC\ , . 

However, the limit yUg ^ of this operator does not exist. This is not accidental; had the 
limit existed, there would be a well-defined operator directly corresponding to the curvature 
F^f^ and we know that even in full LQG, while holonomy operators are well defined, there 
are no operators corresponding to connections and curvatures. This feature is intimately 
intertwined with the quantum nature of Riemannian geometry of LQG. The viewpoint in 
LQC is that the failure of the limit to exist is a reminder that there is an underlying quantum 
geometry where eigenvalues of the area operator are discrete, whence there is a smallest non- 
zero eigenvalue. A, i.e., an area gap [49, 50j. Thus, quantum geometry is telling us that 
it is physically incorrect to let /Iq go to zero because that limit corresponds to shrinking 
the area enclosed by loops Dj^ to zero. Rather, the 'correct' field strength operator in the 
quantum theory should be in fact non-local, given by setting fio in l{2.14\ ) to a non-zero value, 
appropriately related to A. In quantum theory, we simply can not force locality by shrinking 
the loops Djj to zero area. In the classical limit, however, we are led to ignore quantum 
geometry, and recover the usual, local field Fab [iJ, |25| . 

There are two ways of implementing this strategy. In this paper, we will discuss the one 
that has been used in the literature [Ij, llJ, |20| . The second strategy will be discussed in 
J38[ ; it has a more direct motivation and is more satisfactory in semi-classical considerations 
especially in presence of a non-zero cosmological constant. 

In the above discussion, fXo enters through holonomies /i^ ° . Now, it is straightforward 

to verify that (every matrix element of) h)^" is an eigenstate of the area operator Ar = \p\ 
(associated with the face of the elementary cell V orthogonal to the fcth direction): 

Arhi'"\c)= f^4il hi^°\c). (2.18) 



In the first strategy, one fixes /Xo by demanding that this eigenvalue be A = 2\/3n'yipi, the 
area gap [iJ, |50[, so that Ho = 3a/3/2. 
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To summarize, by making use of some key physical features of quantum geometry in 
LQG, we have arrived at a 'quantization' of the classical constraint Cgrav^ 

C'grav = C'(^:i \,._^ (2.19) 



grav 1^^ 

There are still factor ordering ambiguities which will be fixed in the section III B 21 
We will conclude our summary of the quantization strategy by comparing this construc- 
tion with that used in full LQG [ll|3l|5l|. As one would expect, the curvature operator there 
is completely analogous to ()2.14|) . However in the subsequent discussion of the Hamiltonian 
constraint certain differences arise: While the full theory is diffeomorphism invariant, the 
symmetry reduced theory is not because of gauge fixing carried out in the beginning of 
section 111 A 11 (These differences are discussed in detail in sections 4 and 5 of J14].) As a 
result, in the dynamics of the reduced theory, we have to 'parachute' the area gap A from 
the full theory. From physical considerations, this new input seems natural and the strategy 
is clearly viable since that the constraint operator has the correct classical limit. However, 
so far there is no systematic procedure which leads us to the dynamics of the symmetry 
reduced theory from that of the full theory. This is not surprising because the Hamiltonian 
constraint in the full theory has many ambiguities and there isn't a canonical candidate 
that stands out as being the most satisfactory. The viewpoint in LQG is rather that, at 
this stage, it would be more fruitful to study properties of LQG constraints such as the one 
introduced in this section and in |3a], and use the most successful of them as a guide to 
narrow down the choices in the full theory. 

2. Constraint operators and their properties 

It is easy to verify that, although the classical constraint function Cgrav we began with 
is real on the phase space Lgj.^^^, the operator Cgra-v is not self-adjoint on H^^^ . This came 
about just because of the standard factor ordering ambiguities of quantum mechanics and 
there are two natural re-orderings that can rectify this situation. 

First, we can simply take the self-adjoint part of (j2.19|) on H^^^ and use it as the gravi- 
tational part of the constraint: 



C' =i 

grav r) 



'--'grav ~r (^tvgj-avj yA.Av) 



It is convenient to express its action on states \l/(/i) := (^|/i) in the fi — or the 
triad/geometry — representation. The action is given by 

C;,, ^ifi) = /; ^i^^ + 4/.,) + /: vI/(/i) + /: vI/(^ - 4/Xo) (2.21) 

where the coefficients /4, fo are functions of /i: 



f'oif^) 



I Air ipi 



3 (7/^0 



3 3 



fW = -^(/:(/i) + /> + 4/.„)) 

/:(/i) = -^(/:(/i) + /:(/i-4/i„)) (2.22) 
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By construction, this operator is self-adjoint and one can show that it is also bounded above 
(in particular, (\1/, Cg^av ^) < a/vt/S ipi (7;Uo)~ 2 (5 2 — 32)(\1'|\1')). Next, consider the 'parity' 
operator 11 defined by n\l'(/x) := \1'(— /i). It corresponds to the fiip of the orientation of 
the triad and thus represents a large gauge transformation. It will play an important role 
in section ITVl Here we note that the functions /4, fo are such that the constraint operator 
C', commutes with 11: 

grav 

[C^U'n] = (2.23) 

Finally, for ^ ^ ^^ = 3-\/3/2, this difference operator can be approximated, in a well- 
defined sense, by a second order differential operator (a connection-dynamics analog of 
the WDW operator of geometrodynamics). Since C'^^^.^ is just the self-adjoint part of the 
operator used in 1^, we can directly use results of section 4.2 of that paper. Set p = 



STT'jfiipi/G and consider functions \l/(p) which, together with their first four derivatives are 
bounded. On these functions, we have 

6471^ d^ d^ 

c^u np) ^ ^- 4i [v^^ + ^v^] np) =■■ c^;.r^ np) (2.24) 

where ~ stands for equality modulo terms of the order 0{^o)- That is, had we left /Xq as a 
free parameter, the equality would hold in the limit /Xq — > 0. This is the limit in which the 
area gap goes to zero, i.e., quantum geometry effects can be neglected. 

Let us now turn to the second natural factor ordering. The form of the expression ()2.17p 
of Cgrav suggests [5^ that we simply 're-distribute' the sin^ /IqC term in a symmetric fashion. 
Then this factor ordering leads to the following self-adjoint gravitational constraint: 



^ 242 (sgnp) 



^■•"^ 87r73/i3£2 



PI 



sm fioC I sm -— V cos — cos -— V sm -— 1 sm fioC 

\ ^ Zj Zj Zj ' 



Mo— 2 



sin /ioC A sin /ioC (2.25) 

For concreteness in most of this paper we will work with this form of the constraint and, 
for notational simplicity, unless otherwise stated set /Iq = 3-\/3/2. However, numerical 
simulations have been performed also using the constraint C*gi.av of Eq. ()2.2()j) and the 
results are robust. 

Properties of C*grav which we will need in this paper can be summarized as follows. First, 
the eigenbasis |/i) of p diagonalizes the operator A. Therefore in the p-representation, A 
acts simply by multiplication. It is easy to verify that 

ivl>(^) = -2./^-^ |^ + ^^|f _|^_^^|f vl>(^) (2.26) 

V 6 (^^,2)5 

By inspection, A is self-adjoint and negative definite on Ti.^^^ . The form of Cgrav now implies 
that it is also self-adjoint and negative definite. Its action on states ^I/(/i) is given by 

C'grav ^(/i) = /+ ^(/i + 4;Uo) + fo^{l2) + /- ^(/X - 4/X,) (2.27) 
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where the coefficients f±, fo are again functions of /x: 

^ V b {^jj2)2 

/o(/i) = -/+(m) -/-(/")• (2.28) 

It is clear from Eqs ()2.25p and ()2.26|) that Cgrav commutes with the parity operator 11 
which ffips the orientations of triads: 

[C^grav, n]=0. (2.29) 

Finally, the 'continuum limit' of the difference operator Cgrav yields a second order differential 
operator. Let us ffist set 

x(/i) = /+(/x - 2fio)[^{fi + 2/i,) - ^{fi - 2/i„)] (2.30) 

Then, 

C'grav*(/i) = x(/i + 2/io) - x(/i - 2/io) (2.31) 

Therefore, if we again set p = 87r7/i£pi/6 and consider functions \E'(p) which, together with 
their first four derivatives are bounded we have 



1287r2 „^ d d^ _ 

^d^^^ d^^ "■ ^^'^^ 



C^grav nP) ^ -^ 4l [:T-V^-T-] =■■ Cav ^ (p) (2-32) 



where ~ again stands for equality modulo terms of the order 0{fio)- That is, in the limit 
in which the area gap goes to zero, i.e., quantum geometry effects can be neglected, the 
difference operator reduces to a WDW type differential operator. (As one might expect, the 
limiting WDW operator is independent of the Barbero-Immirzi parameter 7.) We will use 
this operator extensively in section IIIII 

In much of computational physics, especially in numerical general relativity, the funda- 
mental objects are differential equations and discrete equations are introduced to approx- 
imate them. In LQC the situation is just the opposite. The physical fundamental object 
is now the discrete equation (J2.27J) with fXo = 3v^/2. The differential equation is the ap- 
proximation. The leading contribution to the difference between the two — i.e., to the error 
— is of the form 0{fio)'^'^ where the second term depends on the wave function \Ef under 
consideration. Therefore, the approximation is not uniform. For semi-classical states, \Ef 
can be large, of the order of 10^^/yU^ in the examples considered in section |Vl In this case, 
the continuum approximation can break down already at /i ~ 10^. 

As one might expect the two differential operators, C*g^v^ and C^^ differ only by a 
factor ordering: 

{c'zr - c^T) m = - ^ 4ibi- m ■ (2.33) 

Remark: As noted in section Hi B 11 the 'continuum limit' fio —^ of any of the quantum 
constraint operators of LQC does not exist on Tif^^^ because of the quantum nature of 
underlying geometry. To take this limit, one has to work in the setting in which the quantum 
geometry effects are neglected, i.e., on the WDW type Hilbert space L^(R, dc). On this 
space, operators Cg^^,^, Cgrav, C^^ and C'^^^^ are all densely defined and the limit can be 
taken on a suitable dense domain. 
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C. Open issues and the model 

In section Hi Al we recalled the kinematical framework used in LQC and in lll BI we extended 
the existing results by analyzing two self-adjoint Hamiltonian constraint operators in some 
detail. Physical states \I'(/i, 0) can now be constructed as solutions of the Hamiltonian 
constraint: 

(C'grav + C'^att)^(/i,0) = O (2.34) 

where (f) stands for matter fields. Given any matter model, one could solve this equation 
numerically. However, generically, the solutions would not be normalizable in the total 
kinematic Hilbert space 'hO^^'^ of gravity plus matter. Therefore, although section IIIBI 
goes beyond the existing literature in LQC, one still can not calculate expectation values, 
fluctuations and probabilities — i.e., extract physics — knowing only these solutions. 
To extract physics, then, we still have to complete the following tasks: 

• In the classical theory, seek a dynamical variable which is monotonically increasing on all 
solutions (or at least 'large' portions of solutions). Attempt to interpret the Hamiltonian 
constraint (j2.34|) as an 'evolution equation' with respect to this 'internal time'. If successful, 
this strategy would provide an 'emergent time' in the background independent quantum 
theory. Although it is possible to extract physics under more general conditions, physical 
interpretations are easier and more direct if one can locate such an emergent time. 

• Introduce an inner product on the space of solutions to ()2.34|) to obtain the physical 
Hilbert space Tiphy The fact that the orientation reversal induced by H is a large gauge 
transformation will have to be handled appropriately. 

• Isolate suitable Dirac observables in the classical theory and represent them by self-adjoint 
operators on Tiphy 

• Use these observables to construct physical states which are semi-classical at 'late times', 
sharply peaked at a point on the classical trajectory representing a large classical universe. 

• 'Evolve' these states using (J2.34J) . Monitor the mean values and fluctuations of Dirac 
observables. Do the mean values follow a classical trajectory? Investigate if there is a drastic 
departure from the classical behavior. If there is, analyze what replaces the big-bang. 

In the rest of the paper, we will carry out these tasks in the case when matter consists of a 
zero rest mass scalar field. In the classical theory, the phase space T^^^ is now 4-dimensional, 
coordinatized by (c, p; (f),P(f,). The basic (non- vanishing) Poisson brackets are given by: 

{c, p} = —^ , and {0, p^} = 1. (2.35) 

The symmetry reduction of the classical Hamiltonian constraint is of the form 

o 

Cg,av + C^matt = ~^ c' ^\ + 8ttG ^ = . (2.36) 

7 Ipp 

Using this constraint, one can solve for c in terms of p and p^. Furthermore, since (p 
does not enter the expression of the constraint, p^ is a constant of motion. Therefore, each 
dynamical trajectory can be specified on the 2-dimensional (p, 0) plane. Typical trajectories 
are shown in Fig. 1. Because the phase space allows triads with both orientations, the 
variable p can take both positive and negative values. At p = the physical volume of the 
universe goes to zero and, if the point lies on any dynamical trajectory, it is an end-point of 
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FIG. 1: Classical phase space trajectories are plotted in the (p,p ^ fj, -plane. For /U > 0, there is 
a branch which starts with a big-bang (at /x = 0) and expands out and a branch which contracts 
into a big crunch (at /x = 0). Their mirror images appear in the /i < half plane. 



that trajectory, depicting a curvature singularity. As the figure shows, for each fixed value 
of p^, there are four types of trajectories, two in the p > half plane and two in the p < 
half plane. Analytically they are given by: 



± 






(2.37) 



where p^, (f)^ are integration constants. These trajectories are related by a 'parity trans- 
formation' on the phase space which simply reverses the orientation of the physical triad. 
As noted before, since the metric and the scalar field are unaffected, it represents a large 
gauge transformation. Therefore, it suffices to focus just on the portion p > of figure 1. 
Then for each fixed value of p^, there are two solutions passing through any given point 
(0^,p^). In one, the universe begins with the big-bang and then expands and in the other 
the universe contracts into a big crunch. Thus in this model, every classical solution meets 
the singularity. 

Finally, we can introduce a natural set of Dirac observables. Since p^ is a constant of 
motion, it is obviously one. To introduce others, we note that is a monotonic function 
on each classical trajectory. Furthermore, in each solution, the space-time metric takes the 
form ds"^ = —dt^ + Vo \p\{t)dSl and the time dependence of the scalar field is given by 



d0 IQiiGpl /— — 

— - = ^ exp ±v 127rG 

dt |p,|f 



(0-0.)] 



(2.38) 



where p^, 0^ and p*. are constants. Thus, in every solution is a monotonic function of time 
and can therefore serve as a good 'internal clock'. This interpretation suggests the existence 
of a natural family of Dirac observables: p\^^, the value of p at the 'instant' 0o. The set 
(p0, p\cj)^) constitutes a complete set of Dirac observables since their specification uniquely 
determines a classical trajectory on the symmetry reduced phase space F*^, i.e., a point in 
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the reduced phase space T^. While the interpretation of as 'internal time' motivates this 
construction and, more generally, makes physics more transparent, it is not essential. One 
can do all of physics on T^: Since physical states are represented by points in T^ and a 
complete set of observables is given by (p^,^!,/,^), one can work just with this structure. 

Remark: In the open i.e., k = model now under consideration, since p is also monotonic 
along any classical trajectory, 0^^ is also a Dirac observable and p could also be used as an 
'internal time'. However, in LQC the expression of the gravitational part of the constraint 
operator makes it difficult to regard p as the 'emergent time' in quantum theory. Moreover 
even in the classical theory, since the universe expands and then recoUapses in the k = 1 
case, p fails to be monotonic along solutions and cannot serve as 'internal time' globally. 

III. WDW THEORY 

The WDW theory of the model has been analyzed in some detail within geometrody- 
namics (see especially J53|). However, that analysis was primarily in the context of a WKB 
approximation. More recently, the group averaging technique was used to construct the 
physical Hilbert space in a general cosmological context J35|, |5J| , an elementary example of 
which is provided by the present model. However, to our knowledge a systematic completion 
of the program outlined in section 111 CI has not appeared in the literature. 

In this section we will construct the WDW type quantum theory in the connection dy- 
namics. This construction will serve two purposes. First it will enable us to introduce the 
key notions required for the completion of the program in a familiar and simpler context. 
Second, we will be able to compare and contrast the results of the WDW theory and LQC 
in detail, thereby bringing out the role played by quantum geometry in quantum dynamics. 

A. Emergent time and the general solution to the WDW equation 

Recall that the phase space of the model is 4-dimensional, coordinatized by (c, p; 0, p^) 
and the fundamental non- vanishing Poisson brackets are given by ()2.35p . The Hamiltonian 
constraint has the form: 

C^r.. + C^ = -^c^V¥\+^^G^ = 0. (3.1) 

7 |p|2 

To make comparison with the standard geometrodynamical WDW theory, it is most con- 
venient to work in the p, representation. Then, the kinematic Hilbert space is given by 
'^kin™ = L^(M^,dpd0). Operators p, operate by multiplication while c and p^ are repre- 
sented as: 

cxiJ = th^—— and p^m = -ih— (3.2) 

Note that, while in geometrodynamics the scale factor is restricted to be non-negative, here 
p ranges over the entire real line, making the specification of the Hilbert space and operators 
easier. 

To write down the quantum constraint operator, we have to make a choice of factor 
ordering. Since our primary motivation behind the introduction of the WDW theory is to 
compare it with LQC, it is most convenient to use the factor ordering that comes from 
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the continuum limit ()2.32|) of the constraint operator of LQC Then, the WDW equation 
becomes: 

-{^-KGhf^^-^ = 8nGh^ \p\-l ^ =: 87rGh^B{p)i^ , (3.3) 



3 



where we have denoted the eigenvalue \p\~^ of |p|^2 by B(p) to facilitate later comparison 
with LQC^ The operator on the left side of this equation is self-adjoint on L^(R, dp) and the 
equation commutes with the orientation reversal operator n\l/(p, 0) = "^{—p, (p) representing 
a large gauge transformation. Thus, if \1/ is a solution to Eq. ()3.3|) . so is "^{—p, (p). 

For a direct comparison with LQC, it is convenient to replace p with /x defined hj p = 
(Sii'yGh/G) /i. Then, ()3.3|) becomes: 

=: -m (3.4) 

where we have recast the equation in such a way that operators involving only (p appear on 
the left side and operators involving only n appear on the right. The WDW equation now 
has the same form as the Klein-Gordon equation in a static space-time, (p playing the role 
of time and B of the (elliptic operator constructed from the norm of the Killing field and 
the) spatial Laplacian. Thus, the form of the quantum Hamiltonian constraint is such that 
(p can be interpreted as emergent time in the quantum theory. In this factor ordering of 
the constraint operator, which emerged in the continuum limit ()2.32j) of LQC, it is not as 
convenient to regard p ~ yU as emergent time. 

Physical states will be suitably regular solutions to ()3.4|) . Since LI is a large gauge 
transformation, we can divide physical states into eigenspaces of P. Physical observables 
will preserve each eigenspace. Since n^ = 1, there are only two eigenspaces, one representing 
the symmetric sector and the other, anti-symmetric. Since the standard WDW theory deals 
with metrics, it is completely insensitive to the orientation of the triad. Therefore, it is 
natural to work with the symmetric sector. Thus, the physical Hilbert space will consist of 
suitably regular solutions '^{fi,(p) to |^^ which are symmetric under fi — * —fi. 

The mathematical similarity with the Klein Gordon equation in static space-times im- 
mediately suggests a strategy to obtain the general solution of ()3.4p . We first note that 
(d/d/i) Y^/I (d/d/i) is a negative definite, self-adjoint operator on L|(M, rf/i), the sym- 
metric sector of L^(M, rf/i). Therefore, G is a positive definite, self-adjoint operator on 
L|(M, B(yu)(i/i). Its eigenf unctions provide us with an orthonormal basis. It is easy to verify 
that the eigenvectors e^(/i) can be labelled by A; G M and are given by 

1 

I 4 



W_^ikln\,i\ 
^k\t^J ■ — 

Their eigenvalues are given by: 



e,(/i) := ^e^^'-^l'^l (3.5) 



ee,(A^) = u;2g^(/i), with ^2 ^ 1^(A;2 + 1) ; (3.6) 



Here, and in what follows, quantities with an underbar will refer to the WDW theory. 
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(where the factor of 1/16 is an artifact of the factor ordering choice which we were led to 
from LQC). The eigenf unctions satisfy the orthonormahty relations: 

d/iB(/i)efc(/i)efc,(/i) = 6{k,k') , (3.7) 

(where the right side is the standard Dirac distribution, not the Kronecker symbol as on 
L^(MBohr5 c^/^Bohr)); and the completeness relation: 

/oo 
dfJ''B{fi)e^{fi) ^(/u) = \fk iff ^(/i) = (3.8) 

-oo 

for any \l/(/i) G L|(R, B(/i)rf/i). 

With these eigenfunctions at hand, we can now write down a 'general' symmetric solu- 
tion to ()3.3p . Any solution, whose initial data at (J) = cpo is such that /i~^/^\l/(/x, 0o) and 
/i~^/''\l/(/i, 0o) are symmetric and lie in the Schwartz space of rapidly decreasing functions, 
has the form: 

/oo 
rfA;§+(%,(/i)e^-<^ + §_(fc)e,(/i)e-^-^ (3.9) 

-oo 

for some \E'±(/i;) in the Schwartz space. Following the terminology used in the Klein-Gordon 
theory, if ^±{k) have support on the negative fc-axis, we will say the solution is 'outgoing' (or 
'expanding') while if it has support on the positive k axis, it is 'incoming' (or, 'contracting'). 
If ^_(/c) vanishes, the solution will be said to be positive frequency and if ^-\-{k) vanishes, 
it will be said to be of negative frequency. Thus, every solution (j3.9|) admits a natural 
decomposition into positive and negative frequency parts. Finally we note that positive 
(respectively negative) frequency solutions satisfy a first order (in (p) equation which can be 
regarded as the square-root of ()3.4|) : 

T^^ = VlvI;^ (3.10) 

where ^/Q is the positive, self-adjoint operator defined via spectral decomposition of O on 
L^(]R, B(/i)(i/i). Regarding as time, this is just a first order Schrodinger equation with a 
non-local Hamiltonian y/Q_. Therefore, a general 'initial datum' f±{n) at = 0o can be 
'evolved' to obtain a solution to ()3.10|) via: 

vl>±(/x,0) = e±^Vi(<^-*°)/±(/x,0„) (3.11) 

B. The physical sector 

Solutions ()3.9jl are not normalizable in TC^^ (because zero is in the continuous part of 
the spectrum of the WDW operator). Our first task is to endow the space of these physical 
states with a Hilbert space structure. There are several possible avenues. We will beg in 
with one that is somewhat heuristic but has direct physical motivation. The idea (55l . l5fil | is 
to introduce operators corresponding to a complete set of Dirac observables and select the 
required inner product by demanding that they be self-adjoint. In the classical theory, such 
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a set is given by p<^ and n\^^. Since p^ commutes with the WDW operator in ()3.4|) . given a 
(symmetric) solution \l/(/i,0) to ()3.4p . 

p<^^(/i,0):=-^;i— (3.12) 

is again a (symmetric) solution. So, we can just retain this definition of p,^ from Ti^in™- The 
Schrodinger type evolutions ()3.11|) enable us to define the other Dirac observable |/i !</,„, where 
the absolute value suffices because the states are symmetric under 11. Given a (symmetric) 
solution \E'(/i, 0) to ()3.4|) . we can first decompose it into positive and negative frequency 
parts \l/±(/i, 0), freeze them at = 0o, multiply this 'initial datum' by |/i| and evolve via 
dSHI): 

i;;:j;„ m{^i, 0) = e^v^(^-*°) |/,| vi/+(/,, 0„) + e-'V0(^-*°) |/,| ii/_(;,, 0„) (3.13) 

The result is again a (symmetric) solution to (|3.4|) . Now, we see that both these operators 
have the further property that they preserve the positive and negative frequency subspaces. 
Since they constitute a complete family of Dirac observables, we have supers election. In 
quantum theory we can restrict ourselves to one superselected sector. In what follows, for 
definiteness we will focus on the positive frequency sector and, from now on, drop the suffix 
+ . 

We now seek an inner product on positive frequency solutions "^{fi, 0) to ()3.4|) (invariant 

under the fi refiection) which makes \p\, and |/i| ,^ self-adjoint. Each of these solutions is 
completely determined by its initial datum \E'(/i, 0o) and the Dirac observables have the 
following action on the datum: 



|/x|0„^(/i,0o) = |/i|^(/i,0o), and ]3<^^(/i,0o) = /iV0^(/^)0o)- (3-14) 

Therefore, it follows that (modulo an overall rescaling,) the unique inner product which will 
make these operators self-adjoint is just: 

(^l|^2)phy=/ dllB{fl)^,^2 (3.15) 

J 4>=4'o 

(see e.g. J55|, |56[). Note that the inner product is conserved, i.e., is independent of the 
choice of the 'instant' = 0o. Thus, the physical Hilbert space 7^"^™ is the space of positive 
frequency wave functions \E'(/i, 0) which are symmetric under fi refiection and have a finite 
norm i3.1^) . The procedure has already provided us with a representation of our complete 
set of Dirac observables on this Ti^^y-. 

i;^„^(/i,0)=e^v^('^-'^°)|^|^(/i,0o), and p<^^(/i, 0) = /iv^^(/i, 0) . (3.16) 



We will now show that the same representation of the algebra of Dirac observables can 
be obtained by the more systematic group averaging method |35l . l54| which also brings out 
the mathematical inputs that go in this choice. (The two methods have been applied and 
compared for a non-trivially constrained system in [57]). Here, one first notes that the total 
constraint operator is self-adjoint on an auxihary Hilbert space Tiau™ ■= -^K^^^! B(/i)d/id0), 
(where, as before the subscript S denotes restriction to functions which are symmetric under 
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/i-reflection) . One must then select an appropriate dense subspace $ of H^^^. A natural 
candidate is the Schwartz space of rapidly decreasing functions. One then 'averages' elements 
of $ under the 1-parameter family group generated by the constraint operator C_ = d'j^ + Q 
to produce a solution to ()H.4|1 : 

7~(fc, u;) e,(/i) e-^ + /(fc, -oo) e,(/i) e— '^) (3.17) 



oo 
oo 



dA; 



oo 



2ku 



where to arrive at the second step we expanded / in the eigenbasis of G and p^ with / as the 
coefficients. Thus, the group averaging procedure reproduces the solution ()3.9|) with ^+{k) = 
f{k, u)/2\uj\ and \E'_(fc) = f{k,—uj)/2\uj\. Solutions \E'/ are regarded as 'distributions' or 
elements of the dual $* of $ and the physical norm is given by the action of this distribution, 
\E'/, on the 'test function' /. However, there is some freedom in the specification of this action 
which generally results in seemingly different but unitarily equivalent representations of the 
algebra of Dirac observables. For us the most convenient choice is: 

/oo /»oo 

d0/ d/iB(/i)'l'^(/i,0)Vl/(/i,0) (3.18) 

-oo J — oo 

where ^E'/(/) is the action of the distribution \1'/ on the test field /. Then, the inner product 
coincides with fl3.15|) and the representation of the Dirac observables is the same as in ()3.12|) 
and ()3.13|) . Had we chosen to drop the factor of a/9 in defining the action of \E'/ on /, we 
would have obtained a unitarily equivalent representation in which the action oi jl\^^ is more 
complicated. 

Finally, with the physical Hilbert space and a complete set of Dirac observables at hand, 
we can now introduce semi-classical states and study their evolution. Let us fix a 'time' 
= 0o and construct a semi-classical state which is peaked at p^ = p*^ and |/i|(^^ = /i*. We 
would like the peak to be at a point that lies on a large classical universe. This implies that 
we should choose jj* ^ 1 and (in the natural classical units c=G=l), p^ ^ h. In the closed 
[k = 1) models for example, the second condition is necessary to ensure that the universe 
expands out to a size much larger than the Planck scale. At 'time' (p = (po, consider the 
state 

dk^{k)ef,{fi)e"^^'l"'-'^*\ where ^(A;) = e — ^^^ . (3.19) 

-oo 

where k* = —y^S/lGnGff p^ and 0* = —^y3/lQnG{\n \fi*\) + (po- It is easy to evaluate the 

integral in the approximation uj = — A/(167rG'/3) k (which is justified because k* <^ —1) 
and calculate mean values of the Dirac observables and their fiuctuations. One finds that, 
as required, the state is sharply peaked at values /i*, p^. The above construction is closely 
related to that of coherent states in non-relativistic quantum mechanics. The main difference 
is that the observables of interest are not /i and its conjugate momentum but rather fi and p^ 
— the momentum conjugate to 'time', i.e., the analog of the Hamiltonian in non-relativistic 
quantum mechanics. 

We can now ask for the evolution of this state. Does it remain peaked at the classical 
trajectory defined by p^ = p^ passing through /i* at = 0o? This question is easy to answer 
because (|3.1ip implies that the (positive frequency) solution \E'(yU, 0) to (|3.4p defined by 
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'initial data' ()3.19|) is obtained simply by replacing 0o by in ()3.19p ! Since the measure of 
dispersion a in ()3.19|) does not depend on 0, it follows that the initial state \l/(/i, 0o) which 
is the semi-classical, representing a large universe at 'time' 0o continues to be peaked at a 
trajectory defined by: 

lnM + 0^. (3.20) 



IQ-kG l/i* 

This is precisely the classical trajectory with p^ = V% passing through /i* at = 0o- This is 
just as one would hope during the epoch in which the universe is large. However, this holds 
also in the Planck regime and, in the backward evolution, the semi-classical state simply 
follows the classical trajectory into the big-bang singularity. (Had we worked the positive 
k* , we would have obtained a contracting solution and then the forward evolution would 
have followed the classical trajectory into the big-crunch singularity.) In this sense, the 
WDW evolution does not resolve the classical singularity. 

Remark: In the above discussion for simplicity we restricted ourselves to eigenfunctions 
§A;(/^) which are symmetric under /z -^ —^ from the beginning. Had we dropped this 
requirement, we would have found that there is a 4- fold (rather than 2-fold) degeneracy 
in the eigenfunctions of 9. Indeed, if 6'(/x) is the step function {9{fj,) = if /i < and 
= 1 if /i > 0), then 9{fi)e\^, ^(/i)e_|^|, 0{—fi)e\^, 6'(— /i)e_|^| are all continuous functions 
of fi which satisfy the eigenvalue equation (in the distributional sense) with eigenvalue 
uj"^ = (167rG/3)(A;^ + 1/16). This fact will be relevant in the next section. 

IV. ANALYTICAL ISSUES IN LOOP QUANTUM COSMOLOGY 

We will now analyze the model within LQG. We will first observe that the form of the 
Hamiltonian constraint is such that the scalar field can again be used as emergent time. 
Since the form of the resulting 'evolution equation' is very similar in the WDW theory, we 
will be able to construct the physical Hilbert space and Dirac observables following the ideas 
introduced in section ITlI Bl 

A. Emergent time and the general solution to the LQC Hamiltonian constraint 

The quantum constraint has the form 

C'grav + C*,^ = (4.1) 

where Cgrav is given by (j2.27|) . Since C*,^ = (SttG) (l/p^/s) (^2-^^ ^^^e constraint becomes: 

SnGpl vl>(/i, 0) = [Bip)]-' C'grav^(/i, 0) (4.2) 

where B{p) is the eigenvalue of the operator l/lp^/^: 



6 \^/^ „_ , „_ / 2 



l6 



^(P)=M^7-^) Bifi), where B{fi)= i—j ^|^ + /i„|4 _ |^ _ ^^jij . (4.3) 

Thus, we now have a separation of variables. Both the classical and the WDW theory 
suggests that could serve as emergent time. To implement this idea, let us introduce 
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an appropriate kinematical Hilbert space for both geometry and the scalar field: "H^"^^ : = 
L^(RBohr, -B(Ai)d/iBohr)®-^^(IR, d0). Smcc is to be thought of as 'time' and /i as the genuine, 
physical degree of freedom which evolves with respect to this 'time', we chose the standard 
Schrodinger representation for but the 'polymer representation' for /i which correctly 
captures the quantum geometry effects. This is a conservative approach in that the results 
will directly reveal the manifestations of quantum geometry; had we chosen a non-standard 
representation for the scalar field, these effects would have been mixed with those arising 
from an unusual representation of 'time evolution'. Comparison with the WDW theory 
would also become more complicated. (However, the use of a 'polymer representation' for (p 
may become necessary to treat inhomogeneities in an adequate fashion.) 
On 7ik°n'^\ the constraint takes the form: 



52^ 



[B{fi)]-' (C+(/x)^(/i + 4/i„, 0) + C°(/i)^(/i, 0) + C-(/i)^(/i - 4/x„, 0)) 
-evl>(/i,0) (4.4) 



were the functions C"^, C° are given by: 



5 



^ (Ai) = ni.. 13 I |y" + 3/ioh - l/U + /io 



3 , ,3 
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TXG 

9|/^o| 
C"(/i) = -C+(/i)-C-(/i). (4.5) 



The form of ()4.4|) is the same as that of the WDW constraint ()3.4|) . the only difference is 
that the 0-independent operator B is now a difference operator rather than a differential 
operator. Thus, the the LQC quantum Hamiltonian constraint can also be regarded as an 
'evolution equation' which evolves the quantum state in the emergent time 0. 

However, since B is a difference operator, an important difference arises from the 
WDW analysis. For, now the space of physical states, i.e. of appropriate solutions to 
the constraint equation, is naturally divided into sectors each of which is preserved by the 
'evolution' and by the action of our Dirac observables. Thus, there is super-selection. Let 
£|e| denote the 'lattice' of points {\e\ + 4n/io, n G Z} on the /i-axis, C-\£\ the 'lattice' of 
points { — \e\ + 4n/io, ra G Z} and let Cs = C\s\ U >C_|e|. Let 'Hf^^'' ,H^^J, and Hf^^ denote 
the subspaces of L^(]RBohr, -S(/i)c?/iBohr) with states whose support is restricted to lattices 
C\s\, jC-ie] and Ce- Each of these three subspaces is mapped to itself by B. Since Cgrav is self- 
adjoint and positive definite on H^^^ = L^(RBohr, d/^Bohr), it follows that B is self-adjoint 
and positive definite on all three Hilbert spaces. 

Note however that since Tif^T^ and H^^^, are mapped to each other by the operator H, 
only Tif^'^ is left invariant by H. Now, because H reverses the triad orientation, it represents 
a large gauge transformation. In gauge theories, we have to restrict ourselves to sectors, 
each consisting of an eigenspace of the group of large gauge transformations. (In QCD in 
particular this leads to the 9 sectors.) The group generated by H is just Z2, whence there 



^ Note that this fundamental evolution equation makes no reference to the Barbero-Immirzi parameter 7. 
If we set /2 = ii/iio and ^{p.,(j)) = 'i{fi,(p), the equation satisfied by ^{jl,(j)) makes no reference to ^io 
either. This is the equation used in numerical simulations. To interpret the results in terms of scale factor, 
however, values of 7 and /Xo become relevant. 
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are only two eigenspaces, with eigenvalues ±1. Since there are no fermions in our theory, 
there are no parity violating processes whence we are led to choose the symmetric sector 
with eigenvalue +1. (Also, in the anti-symmetric sector all states are forced to vanish at 
the 'singularity' /x = while there is no such a priori restriction in the symmetric sector.) 
Thus, we are primarily interested in the symmetric subspace oiTif'^'^; the other two Hilbert 
spaces will be useful only in the intermediate stages of our discussion. 

Our first task is to explore properties of the operator B. Since it is self-adjoint and 
positive definite, its spectrum is non-negative. Therefore, as in the WDW theory, we will 
denote its eigenvalues by u^. Let us first consider a generic e, i.e., not equal to or 2/io. 
Then, on each of the two Hilbert spaces 'H^^J,, we can solve for the eigenvalue equation 
0e^(/i) =u;2e^(/i), i.e., 

C+(/i)e^(/i + 4/i„) + C°(/i)e^(/i) + C-(/i)e^(/i - 4/i„) = u^B{fi)e^{fi) (4.6) 

Since this equation has the form of a recursion relation and since the coefficients C^{fi) never 
vanish on the 'lattices' under consideration, it follows that we will obtain an eigenfunction by 
freely specifying, say, ^(/^*) and \E'(/i* + 4/io) for any fi* on the 'lattice' C\e\ or C-\e\. Hence 
the eigenfunctions are 2-fold degenerate on each of Hf^r and Ti^^T- On Hf^'' therefore, 
the eigenfunctions are 4-fold degenerate as in the WDW theory. Thus, Tif^"^ admits an 
orthonormal basis ef, where the degeneracy index I ranges from 1 to 4, such that 

(e^leS) = Sjj,6{u;,u;'). (4.7) 

(The Hilbert space Hf^^ is separable and the spectrum is equipped with the standard 
topology of the real line. Therefore we have the Dirac distribution 6{uj, u') rather than the 
Kronecker delta S^^^'.) As usual, every element ^(^i) of Tif^^ can be expanded as: 

^(/i)=/" duj^i{uj)e'M where ^/(cu) = (e^l^), (4.8) 

where the integral is over the spectrum of O. The numerical analysis of section |3 and 
comparison with the WDW theory are facilitated by making a convenient choice of this basis 
in Tif^^, i.e., by picking specific vectors from each 4 dimensional eigenspace spanned by e^j. 
To do so, note first that, as one might expect, every eigenvector e^(/i) has the property that it 
approaches unique eigenvectors e'^iuj) of the WDW differential operator 6 as yU — > ±oo. The 
precise rate of approach is discussed in section IXT^ In general, the two WDW eigenfunctions 
e=''(a;) are distinct. Indeed, because of the nature of the WDW operator 6, its eigenvectors 
can be chosen to vanish on the entire negative (or positive) yU-axis; their behavior on the two 
half lines is uncorrelated. (See the remark at the end of section UlIBI ) Eigenvectors of the 
LQC O on the other hand are rigid; their values at any two lattice points determine their 
values on the entire lattice >C±|e|. Second, recall that the spectrum of the WDW operator 9 
is bounded below by u^ > ttG/3, whence e^ with lj"^ < iiG/S does not appear in the spectral 
decomposition of 6. Note however that solutions to the eigenvalue equation Oe^ = i^^e^ 
continue to exist even for u"^ < ttG/3. But such eigenfunctions diverge so fast as /x ^ oo or as 
fi —>■ — oo that (e^|\l/) fails to converge for all "^ e L^(M, B(/i)(i/i), whence they do not belong 
to the basis. What is the situation with eigenvectors of the LQC 9? Since eigenvectors e^ of 
9 approach those of 9, (et^|\&) again fails to converge for all \l/ e TCf^^ if u'^ < nG/S. Thus 
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the spectrum of is again bounded below by ttG/S.^ Therefore, to facihtate comparison 
with the WDW theory we will introduce a variables k via u;^ — vrG/S = {167TG/3)k'^ and use 
k in place of u to label the orthonormal basis. To be specific, let us 

i) Denote by e^|^|(/i) the basis vector in TiS'^T with eigenvalue tu^, which is pro- 
portional to the WDW e_|;i.| as fi ^ oo; (i.e., it has only 'outgoing' or 'expanding' 
component in this limit); 

ii) Denote by em(/i) the basis vector in H^^^, with eigenvalue lu'^ which is or- 
thogonal to e^u,|(/i) (since eigenvectors are 2-fold degenerate in each of 'H^^^,, 
the vector ej^|(/i) is uniquely determined up to a multiplicative phase factor.) 



As we will see in section IV Al this basis is well-suited for numerical analysis. 

We thus have an orthonormal basis e^ in Hf^^ with k G W: (e^|e^,) = 6{k,k'), and 
(^fc kfc') ~ 0- The four eigenvectors with eigenvalue u"^ are now ej^,, e\^ which have support 
on the 'lattice' C\s\, and e|^|,e~|^| which have support on the 'lattice' C^\s\. We will be 
interested only in the symmetric combinations: 

ei'\l^) = I (4(/^) + eti-f^) + e-M + e^(-/i)) (4.9) 

which are invariant under II. Finally we note that any symmetric element ^l'(/i) of Tif^'' can 
be expanded as 

/oo 
dk^>{k)ef{ii) (4.10) 

-oo 

We can now write down the general symmetric solution to the quantum constraint (j4.4p 
with initial data in Tif^'' : 

/oo 
dk [^>^{k)e'i\^^)e'-* + i> .{k)ef {^^)e-^-^] (4.11) 

-oo 

where ^±{k) are in L^(R, dfc). As /i ^ ±cxd, these approach solutions ()H.9|1 to the 
WDW equation. However, the approach is not uniform in the Hilbert space but varies 
from solution to solution. As indicated in section Hi Bl the LQC solutions to ()4.4|) which are 
semi-classical at late times can start departing from the WDW solutions for relatively large 
values of /i, say /i ~ iC^/io- 

As in the WDW theory, if '^ _{k) vanishes, we will say that the solution is of positive 
frequency and if '^+{k) vanishes we will say it is of negative frequency. Thus, every solution 
to (j4.4p admits a natural positive and negative frequency decomposition. The positive 
(respectively negative) frequency solutions satisfy a Schrodinger type first order differential 
equation in 0: 

Ti^ = Ve^± (4.12) 

but with a Hamiltonian v© (which is non-local in fx). Therefore the solutions with initial 
datum ^(yU, 0o) = /±(/^) are given by: 

^±(/i,0) = e±^^('^-*°)/±(/i,0) (4.13) 



^ A rigorous version of this argument can be constructed e.g. by using the Gel'fand triplet [5^ associated 
with the operator Q. However, this step has not been carried out. 
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Remark: In the above discussion, we considered a generic e. We now summarize the 
situation in the special cases, e = and e = 2/io. In these cases, differences arise because 
the individual lattices are invariant under the reffection /i — >■ — /i, i.e., the lattices £|e| and 
£_|e| coincide. As before, there is a 2-fold degeneracy in the eigenvectors of 6 on any one 
lattice. For concreteness, let us label the Hilbert spaces 'Hf^^' and choose the basis vectors 

e~l^{fi), with fc G M as above. Now, symmetrization can be performed on each of these Hilbert 
spaces by itself. So, we have: 

e'i:\^i) = ^{et{f^) + et{-fi)) (4.14) 

However, the vector e;^| (yu) coincides with the vector e_L,{ix) so there is only one symmet- 
ric eigenvector per eigenvalue. This is not surprising: the original degeneracy was 2-fold 
(rather than 4-fold) and so there is one symmetric and one anti-symmetric eigenvector per 
eigenvalue. Nonetheless, it is worth noting that there is a precise sense in which the Hilbert 
space of symmetric states is only 'half as big' in these exceptional cases as they are for a 
generic e. 

For e = 2/io, there is a further subtlety because C"*" vanishes at /i = — 2/io and C~ van- 
ishes at yU = 2/io. Thus, in this case, as in the WDW theory, there is a decoupling and 
the knowledge of the eigenfunction e^(/i) on the positive /z-axis does not suffice to deter- 
mine it on the negative /i axis and vice- versa. However, the degeneracy of the eigenvectors 
does not increase but remains 2-fold because ()4.4j) now introduces two new constraints: 
C±(2/io)e+(±6/io) = [u^B{±2fio) - C°(±2/io)]e+(±2yUo) = 0. Conceptually, this difference 
is not significant; there is again a single symmetric eigenfunction for each eigenvalue. 

B. The Physical sector 

Results of section IIVAI show that while the LQC operator differs from the 
WDW operator in interesting ways, the structural form of the two Hamiltonian constraint 
equations is the same. Therefore, apart from the issue of superselection sectors which arises 
from the fact that is discrete, introduction of the Dirac observables and determination 
of the inner product either by demanding that the Dirac observables be self-adjoint or by 
carrying out group averaging is completely analogous in the two cases. Therefore, we will 
not repeat the discussion of section IIIIBI but only summarize the final structure. 

The sector of physical Hilbert space Hp^y labelled by e G [0, 2/io] consists of positive 
frequency solutions \E'(/i, 0) to ()4.4p with initial data \&(/i, 0o) in the symmetric sector of 
Hgrav- Eq. ()4.11|) implies that they have the explicit expression in terms of our eigenvectors 

^(/i,</))= / dk^{k)e^^\fx)e'^t (4.15) 



where, as before, cj^ = {16nG/3){k'^ + 1/16) and el'\fi) is given by (gS)) and KT^ . By 
choosing appropriate functions "^{k), this expression will be evaluated in section IV Al using 
Fast Fourier Transforms. The resulting \E'(/i, 0) will provide, numerically, quantum states 
which are semi-classical for large /i. The physical inner product is given by: 

(^l|^2). = Yl 5(/i)'l'i(/i,0o)^2(/i,0o) (4.16) 

At6{±|e|+4n/io;neZ} 
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for any 0o. The action of the Dirac observables is independent of e, and has the same form 
as in the WDW theory: 

|/4;vl;(;,,0) = e^v^(*-<^°)|/i|vl>(;.,0„), and p^vl/(;,, 0) = _ ^/^^^^^ . (4.17) 

The kinematical Hilbert space 7^};°*^' is non-separable but, because of super-selection, 
each physical sector Ti-^^y is separable. Eigenvalues of the Dirac observable |/i|^^ constitute 
a discrete subset of the real line in each sector. In the kinematic Hilbert space Tifin^' ^^^ 
spectrum of p is discrete in a subtler sense: while every real value is allowed, the spectrum 
has discrete topology, reflecting the fact that each eigenvector has a finite norm in Ti^in^- 
Thus, the more delicate discreteness of the spectrum of p on H^^^ descends to the standard 
type of discreteness of Dirac observables. Question is often raised whether the kinematic 
discreteness in LQG will have strong imprints in the physical sector or if they will be washed 
away in the passage to the physical Hilbert space. A broad answer — illustrated by the area 
eigenvalues of isolated horizons [59|, |60| — is that the discreteness will generically descend 



to the physical sector at least in cases where one can construct Dirac observables directly 
from the kinematical geometrical operators [l|. The present discussion provides another 
illustration of this situation. 

Note that the eigenvalues of \fi\^^ in distinct sectors are distinct. Therefore which sector 
actually occurs is a question that can be in principle answered experimentally, provided 
one has access to microscopic measurements which can distinguish between values of the 
scale factor which differ by ~ 10£pi. This will not be feasible in the foreseeable future. Of 
greater practical interest are the coarse-grained measurements, where the coarse graining 
occurs at significantly greater scales. For these measurements, different sectors would be 
indistinguishable and one could work with any one. 

The group averaging procedure used in this section is quite general in the sense that it is 
applicable for a large class of systems, including full LQG if, e.g., its dynamics is formulated 
using Thiemann's master constraint program [6l|. In this sense, the physical Hilbert spaces 
Tipj^y constructed here are natural. However, using the special structures available in this 
model, one can also construct an inequivalent representation which is closer to that used in 
the WDW theory. The main results on the bounce also hold in that representation. Although 
that construction appears to have an ad- hoc element, it may well admit extensions and be 
useful in more general models. Therefore, it is presented in Appendix O 

V. NUMERICS IN LOOP QUANTUM COSMOLOGY 

In this section, we will find physical states of LQG and analyze their properties numer- 
ically. This section is divided into three parts. In the first we study eigenfunctions e^(/x) 
of B and then use them to directly evaluate the right side of ()4.15p . thereby obtaining a 
'general' physical state. In the second part we solve the initial value problem starting from 
initial data at = 0^, thereby obtaining a 'general' solution to the difference equation (j4.4j) . 
In the third we summarize the main results and compare the outcome of the two methods. 
Readers who are not interested in the details of simulations can go directly to the third 
subsection. 

A large number of simulations were performed within each of the approaches by varying 
the parameters in the initial data and working with different lattices Cg. They show that 
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the final results are robust. To avoid the making the paper excessively long, we will only 
show illustrative plots. 

A. Direct evaluation of the integral representation ()4.15j) of solutions 

The goal of this sub-section is to evaluate the right side of ()4.15|) using suitable momen- 
tum profiles ^{k). This calculation requires the knowledge of eigenfunctions e^* (/i) of 9. 
Therefore, we will first have to make a somewhat long detour to numerically calculate the 
basis functions e^ (/i) and e^'^ (/x) introduced in section II V Al The integral in (j4.15|) will be 
then evaluated using a fast Fourier transform. 



1. General eigenfunctions of the Q operator: asymptotics 

We will first establish properties of the general eigenfunctions e^uifi) of G that were used 
in section HV Al 

Let us fix a 'lattice', say C\s\. Since the left side of (j4.6j) approaches C^^^e^{fi) as 
|/i| — >• oo, in this limit one would expect each e^(/i) to converge to an eigenfunction e^(/i) 
of G with the same eigenvalue. Numerical simulations have shown that this expectation is 
correct and have also provided the rate of approach. 
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FIG. 2: Crosses denote the values of an eigenfunction et^(^) of for e = and u) = 20. The solid 
curve is the eigenfunction e^(//) of the WDW G to which e^(;u) approaches at large positive fi. As 
/i increases, the set of points on C^ becomes denser and fill the solid curve. For visual clarity only 
some of these points are shown for /x > 100. 

Recall that each WDW eigenfunction e^{fi) is a linear combination of basis functions 
6|fc|(Ai)) §-|fc|(/^) defined in section ITllAI Therefore, given an e^(/i) it suffices to calculate the 
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coefficients of the decomposition of e^{fJ.) with respect to this basis. The method of finding 
these coefficients is presented in detail in Appendix El ^ Once the hmiting e^(;u) were found, 
they were compared with the original eigenfunction ea;(/x) for a variety of values of u. An 
illustrative plot comparing e^(/i) with its limit is shown in Fig. |21 In general, each e^j^fj.) 
approaches distinct eigenf unctions e^ _|_(/i) of 9 in the limits /i -^ ±00. The rate of approach 
is given by: 

_i ,, L-3e^, + (/i) + 0(^) , for/i>0, 

[l^'~'^,-if^) + 0[^^) , for/i<0. 

Numerical tests were performed up to \fi\ = lO^fio- The quantity yU^/'' I ^o; — e^ ± I (/^) was 
found to be bounded. The bound decreases with fi. For the case e = 0, cj = 20 depicted in 
Fig. 121 the absolute bound in the /i-interval (lO^yUo, lO^/Xo) was less than 90. 

Because the eigenf unctions e^(/i) of 9 are determined on the entire lattice C\s\ by their 
values on (at most) two points, the WDW limits for positive and negative fi are not inde- 
pendent. Thus, if the limits are expressed as 

ec^(/i) ^^ Ae|^|(/i) + 5e_|fc|(/i), e^{fi) ^^^^ C e^^^ifi) + D e_\^^{fi) (^21) 

the coefficients C,D are uniquely determined by values of A,B (and vice versa). One 
relation, suggested by analytical considerations involving the physical inner product, was 
verified in detail numerically: 

\A\^-\B\^ = |C|2-|D|2 . (5.3) 

It will be useful in the analysis of basis functions in the next two sub-sections. 

2. Construction of the basis e_,u 

In section HVAI we introduced a specific basis of Tie which is well adapted for comparison 
with the WDW theory. We will now use numerical methods to construct this basis and 
analyze its properties. Our investigation will be restricted to the vectors e^,^{fx) because, 
the physical states \E'(/i,0) of Eq. ()4.15|) we are interested will have negligible projections 
on the vectors ej^|(/i). (Recall that in general \E' G He has support on C\^\ U jC-\£\. e^|^|(/x) 
has support on >C+|e| and e~|^|(/i) on >C_|e|.) 

Each of the eigenfunctions e^,^, is calculated as follows. To solve ()4.6|) . we need to specify 
initial conditions at two points on each of the two lattices. We fix large positive /i^ G C±\s\ 
and demand that the values of e^|^|(/x) agree with those of e_|^|(/i) at the points /i^ and 

/i^ + 4/io. Then e^,^, are evaluated separately on finite domains C±\si fl [— yU±, /i±] of each of 
the two lattices. For large negative /i, these eigenfunctions are linear combinations of the 
WDW basis functions C^e|jr.| + D"^ e_\^. The coefficients C^,D^ are evaluated using the 
method specified in Appendix JBl 



^ If w^ < 7rG/3, then w^ is not part of the spectrum of the self-adjoint operator Q. Nonetheless, by directly 
solving the eigenvalue equation Oc^ — uj'^e^ one can introduce an analogous decomposition onto fixed 
eigenfunctions, e_|_|j,/| := |/^|^^'^ , k' ~ 1/16 — 3/(167rG)(jj^ and also write the limit of e^^ in terms of 
coefficients in this 'basis'. 
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FIG. 3: The exponential growth of |e_|^||(//) in the 'genuinely quantum region' is shown for three 
different values of uj, where uj is given by G e_ 
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;enfunctions e^|^|(/i) were calculated for approximately 2 x 10^ different |/c|'s in the 



Eij 
range 5 < cu < 10"^. They revealed the following properties: 



1 



Each e^|^|(/i) is well approximated by a WDW eigenfunction until one reaches the 
'genuinely quantum region'. In this region the absolute value |e_|, J grows very quickly 



(ii) 



|fc|i 

as n decreases: le^i^il oc e^'^^'^'-^-'^v'^''^', where a = a{e) is a constant on any given 
lattice. This property is illustrated by Fig. |S1 This region of rapid growth is symmetric 
about /i = and its size depends linearly on u (the square-root of the eigenvalue of 
6); its boundary lies at /i ~ 0.5u; /io). (However, this region excludes the interval 
[— 4/io, 4/io] where -B(/i) decreases and goes to zero, departing significantly from its 
WDW analog B(/i).) 

After leaving this region of growth, the basis function e^,^{^) again approaches some 
WDW eigenfunction 

e±^|il51zi,C7±e|,| + Z}±e_|,|. (5.4) 

where the coefficients C^, D^ are large. Their absolute values grow exponentially 
with UJ. To investigate this property qualitatively we defined an 'amplification factor' 



A^ 



\C^\ + \D^ 



(5.5) 



The numerical calculations show that A+ = A , so parts of e^,^, supported on >C+|e| 
and >C-|e| are amplified equally. The dependence of A^ on u and e is shown in Fig. HI 
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FIG. 4: The amplification factor A^ in the 'genuinely quantum region is shown as a function of 
the parameter e labeling the lattice and uj. 
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FIG. 5: The function a(e) of Eo 15. 61 is plotted by connecting numerically calculated data points. 
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FIG. 6: The function b{e) of Eg 15.61 is plotted by connecting numerically calculated data points. 



Almost everywhere it can be well approximated by the function 



X^{u;,e) 



au>+b 



(5.6) 



where a,b are rather complicated functions of e. The fits of a{e),b{e) are presented 
in Fig. |31 and Fig. IHl The actual simulations were carried out for various values of 
p^ = hu, up to p^ = 10^. 

(iii) The general relation ()5.H|1 holds in our case with \B\ = 0. Existence of the tremendous 
amplification now implies that the absolute values of C and D are almost equal (with 
differences of the order of 1/A^) 

\C^\ ^ \D^\ . (5.7) 

Thus for negative /i, eigenf unctions asymptotically approach WDW eigenfunctions 
and are almost equally composed of incoming and outgoing waves. 



3. Basis for the symmetric sector 



Once the basis functions e_|^| are known one can readily use ()4.9|) to construct the basis 

hi 

be 
The numerical calculations show the following properties: 

(i) Each symmetric basis eigenfunction is strongly suppressed in the 'genuinely quantum 
region' around to /i = 0. The behavior of e,,, in this region implies that |e_|,|| ex 



e_L\ for solutions to ()4.6p which are symmetric under /i -^ — yU. Due to strong amplification 

(s) + 

in the region around /i = the behavior of eL, is dominated by properties of e_,u for /i < 0. 



cosh.[a -yyuJI), where a is a function of e only. 



|fc| 



-|fc|i 
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(ii) Outside the 'genuinely quantum region' e_L, quickly approaches the 
WDW eigenfunction almost equally composed of incoming and outgoing 'plane 
waves' (e|;.| and e_|^|). In the exceptional cases e = or 2/io the two contributions are 
exactly equal. To establish this result, note first that the symmetry requirement and 
the fact that C+(— 2/Xo) = C~(2/io) = imply that in both cases the value of e_!u at 
one point already determines the complete 'initial data' (i.e., values at some /x* and 
/i* + 4/io) and hence the eigenfunction on the entire C^'- 



'-|fc| 
'-|fe| 



(6/io) 



M 



Hfcr 



-4Aio) 



'-|fc| 



uj^B{2i2o) - C°(2/i,) 
C+{2fio) 



(0), 

^-\k\ 



C^l^o) , 



for e 
for e 



0, 
2/io • 



(5.8a) 
(5.8b) 



Because of the reality of coefficients of Eq. ()4.fi|l . this implies that the phase of e_L, 
is exactly constant, whence contributions of ei^i and e_|;.| are also exactly equal. 

(iii) On each lattice C±\e\ the incoming and outgoing components are rotated with respect 



to each other by an angle a- 



.W 



/i^OO 



Hfcl 



'±|e 



^^^(e'" §1^1 +e-- e 



&-\k\) 5 



(5.9) 



where z^ are some complex constants satisfying l^;"*"! ?^ |z |, while the phases a^ are 
functions of e and u. In general, for £ 7^ or e 7^ 2fio, a+ need not equal a_. 




1.2*10' 



FIG. 7: Plot of the wave function ^(/U, 0) obtained by directly evaluating the right side of 1)4. 15() . 
Parameters are p^ = 500, Ap^/p^ = 0.05 and e = Ho- 
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FIG. 8: Expectation values and dispersion of \fi\^ for wave function presented in Fig. [7| are 
compared with classical trajectories. 



4- Evaluation of the integral in i4-15\ ) 



.(*) 



Now that we have the symmetric basis functions €,_L, at our disposal, we can obtain the 
desired physical states by directly evaluating the integral in ()4.15j) . 

We wish to construct physical states which are sharply peaked at a phase space point 
on a classical trajectory of the expanding universe at a late time (e.g., 'now'). The form of 
the integrand in (|4.15|) . the expression (|4.16|) of the physical inner product, the functional 
form of a;(/c) and standard facts about coherent andsqueezed states in quantum mechanics 
provide a natural strategy to select an appropriate \1/(A;). If we set 



^{k) 






-iu)4>* 



(5.10) 



with suitably small a, the final state will be sharply peaked at 



levrG/l^ 



(5.11) 



and the parameter 0* will determine the value /i* of the Dirac observable l/^l^^ at which 
the state will be peaked at 'time' (po- As mentioned in section 1111 Bl to obtain a state 
which is semi-classical at a late time, we need a large value of p<^: p^ ^ /i in the classical 

units, c=G=l. Therefore, we need k* <C —1, whence the functions '^{k) of interest will be 
negligibly small for /c > 0. Therefore, without loss of physical content, we can set them 
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to zero on the positive k axis. This is why the exphcit form of eigenfunctions e|^| is not 
required in our analysis. 

Thus, the integral we wish to evaluate is: 

^(/i,0) = / dke—^^el^\fi)e''^^^^^^-^*\ (5.12) 

where a is the spread of the Gaussian. The details of the numerical evaluation can be 
summarized as follows. 

(s) 

• For a generic e, the e_!u were found numerically following the procedure specified in 
sections fV A 21 and IV A 31 For the exceptional cases, e = 0, 2/io, in order to avoid loss 

(s) (s) 

of precision in the region where e_L, is very small, we provided 'initial values' of e^; 
at /i = ±e and /i = ±e + 4/io using ()5.8|) . On the k axis we chose a set {kj} of points 
which are uniformly distributed across the interval [k* — 10a, k* + lOcr]. In numerical 
simulations, the number / of points in the set {kj} ranged between 2^^ and 2^^. 

• Next, for each kj, we calculated e_L,{fii) for /i G {±e + Anfio '■ n G {— A^, . . . , A^}} 
where A^ is a large constant (~ 50pp. 

• Finally, we evaluated ()5.12|) using fast Fourier transform. The result was a set of 
profiles \E'(/ij, 0j) where (j)j = ^J^i-n /{AG){j — l/2)/{ki — ki) This is a positive frequency 
solution to the LQC equation ()4.4j) . 

Our next task is to analyze properties of these solution. Given any one \l/(0j,/ij), we 
chose 'instants of time' and calculated the norm and the expectation values of our Dirac 
observable l/tL and pA, using: 



^ Y, B{^^,)\^^,\\^{(t),^^,)W (5.13b) 

(S>4>) = T^2 H Bifi,)^{<p,fi,) ^d^^{(f),fii). (5.13c) 






Finally, the dispersions were evaluated using their definitions: 



2 / , / V 2 



These calculations were performed for 16 different choices of e and for ten values of p^ 
up to a maximum of Prp = 10^, and for 5 different choices of the dispersion parameter a. 
Fig. [7| provides an example of a state constructed via this method. The expectation values 

of ( |yu|0 ) are shown in Fig. |H1 

Results are discussed in section EO below. 
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B. Evolution in (j) 

We can also regard the quantum constraint ()4.4p as an initial value problem in 'time' 
and solve it by carrying out the 0-evolution. Conceptually, this approach is simpler since 
it does not depend on the properties of the eigenfunctions of G. However, compared to the 
direct evaluation of the integral (j4.15|) . this method is technically more difficult because it 
entails solving a large number of coupled differential equations. Nonetheless, to demonstrate 
the robustness of results, we carried out the 0-evolution as well. This sub-section summarizes 
the procedure. 

1. Method of integration 

At large |/i| the difference equation is well approximated by the WDW equation which 
is a hyperbolic partial differential equation (PDE). However since G couples /i in discrete 
steps, for a given e-sector fl4.4|) is just a system of a countable number of coupled ordinary 
differential equations (DDEs) of the 2nd order. ^ 

Technical limitations require restriction of the domain of integration to a set |/i — el < 
4:Nfio, where 1 -C iV e Z. This restriction makes the number of equations finite. However 
one now needs to introduce appropriate boundary conditions. The fundamental equation on 
the boundary is id^"^ = svG^, where s = +1 (—1) for the forward (backward) evolution in 
(j). It is difficult to calculate vQ^ at each time step. Therefore, just on the boundary itself, 
this equation was simplified to: 

a<^^(/i,0) = sv/^^(|^|-2/x„)(^(/i,0)-^(/x-4sgn(/iK,0)) . (5.15) 

This is the discrete approximation of the continuum operator d^'^ = s{^ylQ^TG/3)^^^'^ 
which itself is an excellent approximation to the fundamental equation when the boundary 
is far. The boundary condition requires the solution to leave the domain of integration. 
(For, to make the evolution deterministic in the domain of interest, it is important to avoid 
waves entering the integration domain from the boundary.) The boundary was chosen to lie 
sufficiently far from the location (in /i) of the peak of the initial wave packet. Its position 
was determined by requiring that the value of the wave function at the boundary be less 
than 10~" times that at its value of the peak, and n ranged between 9 and 24 in different 
numerical simulations. 

Three different methods were used to specify the initial data \l/ and 9<^\l/ at (j) = (po- These 
are described in section IVB21 The data were then evolved using the fourth order adaptive 
Runge-Kutta method (RK4). To estimate the numerical error due to discretization of time 
evolution, two sup-norms were used: 

sup^,^^^^^^_*dW ^ (5,16) 

SUP|^,_e|<7V/.J^2|(0) 



Unfortunately for the £ = sector the equation is singular at /i = 0, so the analysis of this sub-section 
will not go through. This sector was handled by the direct evaluation of the integral representation of the 
solution, presented in the last sub-section. 
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and 



1^1 -*2|// 



(</>) 



SUP|^,_e|<iV^J|^l| - 1^211 (0) 



(5.17) 



Fig. 121 shows an example of the results of convergence tests for the solution corresponding 
to ptf) = 10^ with initial spread in p^ of 3% (used in figures fT^ — fTS|) . One can see that the 
phase of \E' is more sensitive to numerical errors than its absolute value. Therefore, although 
the accuracy is high for the solution \E'(yU, (p) itself, it is even higher for the mean values and 

dispersions of \fi\(f,„ and p^. 




FIG. 9: Error functions |^(Af' )— ^(Af)|/ (upper curve) and |^(m')~^(M)|// (lower curve) are plotted 
as a function of time steps. Here ^im') refers to final profile of wave function for simulation with 
M' time steps, "^(m) refers to the final profile for the finest evolution (approximately 2.88 x 10^ 
time steps). In both cases, the evolution started at </> = and the final profile refers to <j) = —1.8. 



2. Initial data 

Although we are interested in positive frequency solutions, to avoid having to take square- 
roots of 9 at each 'time step', the second order evolution equation (j4.4|) was used. Thus, 
the initial data, consist of the pair \E' and its time derivative d^^^, specified at some 'time' 
(f)o- The positive frequency condition was incorporated by specifying, in the initial data, the 
time derivative of \Ef in terms of "$. Since ^I/(/i) = \1'(— /x), we will restrict ourselves to the 
positive /i axis. The idea is to choose semi-classical initial data peaked at a point (p^, /i*) 

on a classical trajectory at 'time' = 0o, with (p* ^ VGh?, ji* ^ 1), and evolve them. 
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To avoid philosophical prejudices on what the state should do at or near the big bang, we 
specify the data on the expanding classical trajectory at 'late time' (i.e., 'now') and ask it 
to be semi-classical like the observed universe. 

The idea that the data be semi-classical was incorporated in three related but distinct 

ways. 

1. Method I: This procedure mimics standard quantum mechanics. Since (c, /i) are canon- 
ically conjugate on the phase space, we chose \I/(/i) to be a Gaussian (with respect 
to the measure defined by the WDW inner product) and peaked at large /i* and the 
value c* of c at the point on the classical trajectory determined by (/i*,p^, 4>o))'- 

^f \^^:= N \fi\-^ e — 5i^e-*^V^. (5.18) 

where A^ is a normalization constant. The initial value of (90\&|(^^ was calculated using 
the classical Hamilton's equations of motion 



9,*U. = I -,,n(k')f-^^-^^ + ,.'^\ * I,. . (5,19) 

where u* = p*^/h. 

2. Method II: This procedure takes advantage of the fact that ()H.19|1 provides a 
WDW physical state which is semi-classical at late times. The idea is to calculate 
\1/ and d^'^ aX cj) = (f)o and use their restrictions to lattices Ce as the initial data for 
LQC, setting uj = a/IOttG/S \k\ as in the discussion following ()3.19|) . Thus the initial 
data used in the simulations were of the form: 

1 

^U^ = ("AVe-^^^'^e^^*^'^^, (5.20) 



a^vl>|^^ = |_sgn(r)^/i|^a2ln^ + z^(r)lv^U„ (5.21) 

This choice is best suited to comparing the LQC results with those of the WDW theory. 

The spreads a of Method I and a of Method II are related to the initial spread A/i|^^ 
as follows 



a 



^ = V2Afi\^^ . (5.22) 

a 

3. Method III: To facilitate comparison with the direct evaluation of the integral solution 
described in section IV A| a variation was made on Method II. Specifically, in the ex- 
pression ()3.19j) of \1/, the WDW basis eigenfunctions e^(/i) were rotated by multiplying 
them with a k and e dependent phase factor defined in Eq. ()5.9|) 

e_|.| ^ e-^'^%_|,| . (5.23) 

These phases were first found numerically using the method specified in appendix iBl 
and then functions of the form 

a+ = A \n{Bk + C)k + D , (5.24) 
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(where A, B, C, D are real constants) were fitted to the results. After subtraction of 
the 0th and the 1st order terms in the expansion in k around k* (which respectively 
correspond to a constant phase and a shift of origin of 0) the resulting function was 
used to rotate the basis e^(/i) appearing in (j3.19p ) via (j5.23p ). The expression on the 
right side of Eq ()3.19j) and its 0-derivative were then integrated numerically at = (po- 

In the simulations, 15 different values of p^ were used ranging between 10^ and 10^. fi* 
was always greater than 2.5p(^. The dispersion a was allowed to have five different values. 
These simulations involved four different, randomly chosen values of e. 



m 




4.0^ 



FIG. 10: The absolute value of the wave function obtained by evolving an initial data of Method 
II. For clarity of visualization, only the values of |^| greater than 10~^ are shown. Being a physical 
state, ^ is symmetric under // -^ —fi. In this simulation, the parameters were: e = 2/io, p^ = W^, 
and Ap^/p^ = 7.5 x 10~^. 



C. Results and Comparisons 

We now summarize the results obtained by using the two constructions specified in sec- 
tions |^^ and |^^ The qualitative results are robust and differences lie in the finer struc- 
ture. In particular, evolutions of initial data constructed from the three methods described 
in section IV B 21 vield physical states \l/(/i,0) which are virtually identical except for small 
differences in the behavior of relative dispersions of the Dirac observables. The numerical 
evaluation of the integral ()4.15|) in section IV A 41 yielded results very similar to those ob- 
tained in IVB2I using initial data of Method III. An example of results is presented in Fig. 
CniandFig. HU 

Highlights of the results can be summarized as follows: 
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FIG. 11: The expectation values (and dispersions) of |/i|(^ are plotted for the wavefunction in Fig 
^land compared with expanding and contracting classical trajectories. 



(1) The state remains sharply peaked throughout the evolution. However, as shown in 
Fig. E[ while the product A0 Ap^ is nearly constant for large fi, there is a substantial 
increase near /i = 0. 

(2) The expectation values of fifp and p^ are in good agreement with the classical trajec- 
tories, until the increasing matter density approaches a critical value. Then, the state 
bounces from the expanding branch to a contracting branch with the same value of 
{Prf,)- (See Fig. ITT|) . This phenomena occurs universally, i.e., in every e-sector, for 
all three methods of choosing the initial data and for any choice of p^ ^ vGh?. In 
this sense the classical big-bang is replaced by a quantum bounce. Note that this is in 
striking contrast with the situation with the WDW theory we encountered in section 
IIIIBl even when the initial data is chosen using Method II which is tailored to the 
WDW theory. As indicated in Appendix El the existence of the bounce can be heuris- 
tically understood from an 'effective theory'. The detailed numerical work supports 
that description, thereby providing a justification for the approximation involved. 



(3) If the state is peaked on the expanding branch /i(0) = /i* exp(^167rG'/3 (0 — 0o)) 
in the distant future, due to the bounce it is peaked on a contracting branch in the 
distant past, given by yu(0) = D{p^) /i*exp(— -^167rG/3 {(f) — (t)o))-i where D{p^) = 
filp'^/12'7iGh'^. Thus, for large \fi\ the solution \&(/x, 0) exhibits reflection symmetry 
(about = 0o — I ^^D{Pa))- However, it is not exactly reflection symmetric (compare 

b])- 
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FIG. 12: Comparisons between the relative dispersions A/i//i as functions of ( 
specifying initial data and the result of direct construction. 



for all three methods 



(4) As a consistency check, we verified that the norm and the expectation value {p^) are 
preserved during the entire evolution. Furthermore, the dispersion also remains small 
throughout the evolution, although the precise behavior depends on the method of 
specification of the initial data. Differences arise primarily near the bounce point and 
manifest themselves through the behavior of the relative dispersion A/i//i as a function 
of 0. These are illustrated (for all three methods as well as for the direct evaluation of 
section IV A|) in Fig. ^1 Finally, as argued in Appendix ^ since the state is sharply 
peaked, the value of Ayu//i can be related to A0 via Eq IA13I One finds that the 
product A0Ap<^ is essentially constant in the region away from the bounce but grows 
significantly near the bounce. 

The differences can be summarize as follows: 

(i) In Method I the initial state is a minimum uncertainty state in (/x, c) but doesn't 
minimize the uncertainty in {(p^Ptp) at any value of 0. The relative spread in ix 
remains approximately constant in the regions where {(x^) is large and increases 
near the bounce point monotonically. The wave function interpolates 'smoothly' 
between expanding and contracting branches (see Fig. IT^ . On the other hand 
the product A(j)Ap^ of uncertainties has a value much higher than 1/2, grows 
quickly near the bounce and settles down to constant value after it. See Fig. UHl 

(ii) The state obtained by evolving the initial data constructed from Method II has 
minimal uncertainty in {(j),p^). A/i//i is approximately constant for large (/IJ) 
and it decreases quickly near the bounce point, reaching its minimal value shortly 
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FIG. 13: A zoom on the absolute value of the wave function near the bounce point. Initial data 
was specified using Method I. 



0.5 




FIG. 14: A zoom on the absolute value of the wave function near the bounce point. Initial data 
was specified using Method II. 



before the bounce point. After the bounce it grows and stabilizes at the value 
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FIG. 15: A zoom on the absolute value of the wave function near the bounce point. Initial data 
was specified using Method III. 
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FIG. 16: The uncertainty product A</>Ap0 for three methods of specifying the initial data. 



of the relative spread found for data constructed using Method I (for the same 
values of p^ and initial A/i//i). Behavior of the wave function is also different 
from that in the previous case. Near the bounce point its value grows to form a 



42 



bulge (see Fig. IT^ . The product A(j)Ap^ remains almost constant for as long as 
the results of the numerical measurement are reliable (see Fig. ITBj) approaching a 
constant after the bounce, with a somewhat higher value. The heuristic estimate 
on AcpApfj, (see Appendix |X)) agrees with these results. Finally, we also found 
that the increase of the relative spread depends on p^ and initial Afi/fi. 

(iii) The state obtained by evolving the initial data constructed from Method III does 
not have minimum uncertainty in {4>,P(j)). The behavior of Afi/fi is similar to the 
previous case except that it becomes equal to Ap^/p^ at the bounce point and its 
asymptotic value in the contracting branch is same as its starting value. Thus the 
spread is symmetric with respect to reflection in around the bounce point. The 
difference between the value of A/i//i for large (p^) and the one corresponding to 
the minimum uncertainty state (for the same p^ and Ap^) is a function of p^ and 
Ap^. The wave function forms a symmetric bulge near the bounce point and the 
value A0Ap0 remains constant within regime of validity of its estimation (see 
FigHni). 

(iv) The direct construction of section IV Al vields results similar to those obtained by 
using the third method to choose the initial data for the evolution. 

(v) The differences between the relative dispersion Afi/fi resulting from different 
methods of choosing initial data can also be estimated using the solutions to the 
effective dynamical equations (for details of the method see Appendix IX|) . 

We will conclude with two remarks. 

(i) Let us return to the comparison between the results of LQC and the WDW theory in 
light of our numerical results. As remarked in section Hi Bl on functions "^{fi) which 
we used to construct the semi-classical initial data, the leading term in the difference 
(C'grav — C'l^av )^ bctwcen the actions of the two constraint operators goes as Oi^^l)"^"". 
Now, in LQC fXo is fixed, (jJo = 3-\/3/2), and on semi-classical states \1/"" ~ /c*^//x'^\l/. 
Hence the difference is negligible only in the regime fc*^//i^ -C 1. In our simulations, 
k* ~ 10^ whence the differences are guaranteed to be negligible only for /i ^ 10^, 
i.e., well away from the bounce. But let us probe the situation in greater detail. Let 
us regard fio as a mathematical parameter which can be varied and shrink it. In the 
limit fio -^ 0, (Cgrav — C'l^av )^ should tend to zero. Since there is no bounce in the 
WDW theory, we are led to ask: Would the LQC bounce continue to exist all the way 
to yUo = or is there a critical value at which the bounce stops? The answer is that for 
any finite value of fio there is a bounce. However, if we keep the physical initial data 
the same, we find that as we decrease /Xq the solution follows the classical trajectory 
into the past more and more and bounce is pushed further and further in to the past. 
In the limit as fio goes to zero, the wave function follows the classical trajectory into 
infinite past, i.e., the bounce never occurs. This is the sense in which the WDW result 
is recovered in the limit /lo ^ 0. 

(ii) Numerical simulations show that the matter density at the bounce points is inversely 
proportional to the expectation value {p^) = p*^ of the Dirac observable pf Given two 
semi-classical states with {p^) = pj and p^, we have Pcrit/Pcrit = P^/p},- Therefore, 
this density can be made small by choosing sufficiently large p^. Physically, this is 
unreasonable because one would not expect departures from the classical theory until 
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matter density becomes comparable to the Planck density. This is a serious weakness 
of our framework. Essentially every investigations within LQC we are aware of has this 
— or a similar — drawback but it did not become manifest before because the physics 
of the singularity resolution had not been analyzed systematically. The origin of this 
weakness can be traced back to details of the construction of the Hamiltonian constraint 
operator, specifically the precise manner in which the operator corresponding to the 
classical field strength F^^ was introduced. The physical idea that in LQC the operator 
corresponding to the field strength F^^, should be defined through holonomies, and that 
quantum geometry does not allow us to shrink the loop to zero size, seem compelling. 
However, the precise manner in which the value of fio was determined using the area 
gap A is not as systematic and represents only a 'first stab' at the problem. In J3a] we 
will discuss an alternate and more natural way of implementing this idea. The resulting 
Hamiltonian constraint has a similar form but also important differences. Because of 
similarities the qualitative conclusions of this analysis — including the occurrence of 
the quantum bounce — are retained but the differences are sufficiently important to 
replace the expression of the critical density by p(.j,j^ = (3/87rG'7^A) where, as before, 
A = 2-\/37r7£pi is the area gap. Since p'^^^^ is of Planck scale and is independent of 
parameters associated with the semi-classical state, such as p^, the departures from 
classical theory now appear only in the Planck regime. This issue is discussed in detail 
in ET 



VI. DISCUSSION 

We will first present a brief summary and then compare our analysis with similar con- 
structions and results which have appeared in the literature. However, since this literature 
is vast, to keep the discussion to a manageable size, these comparisons will be illustrative 
rather than exhaustive. 

In general relativity, gravity is encoded in space-time geometry. A basic premise of LQG 
is that geometry is fundamentally quantum mechanical and its quantum aspects are central 
to the understanding of the physics of the Planck regime. In the last three sections, we 
saw that LQC provides a concrete realization of this paradigm. In our model every classical 
solution is singular and the singularity persists in the WDW theory. The situation is quite 
different in LQC. As in full LQG, the kinematical framework of LQC forces us to define 
curvature in terms of holonomies around closed loops. The underlying quantum geometry 
of the full theory suggests that it is physically incorrect to shrink the loops to zero size 
because of the 'area gap' A. This leads to a replacement of the WDW differential equation 
by a difference equation whose step size are dictated by A. Careful numerical simulations 
demonstrated in a robust fashion that the classical big-bang is now replaced by a quantum 
bounce. Thus, with hindsight, one can say that although the WDW theory is quite similar 
to LQC in its structure, the singularity persists in the WDW theory because it ignores the 
quantum nature of geometry.^ 



^ Differences arise in two places. The first occurs in the matter part of the constraint and stems from the 
fact that the functions B{^) representing the eigenvalues of the operator l/jpp/^ in LQG is different from 
the corresponding B(/i) of the WDW theory. The second comes from the role of quantum geometry in 
the gravitational part of the Hamiltonian constraint, emphasized above. In our model, qualitatively new 
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Emergent time 

In this paper we isolated the scalar field (f) as the emergent time and used it to motivate 
and simplify various constructions. However, we would like to emphasize that this — or 
indeed any other — choice of emergent time is not essential to the final results. In the 
classical theory (for any given value of the constant of motion p^), we can draw dynamical 
trajectories in the /^ — </> plane without singling out an internal clock. A complete set of 
Dirac observables can be taken to be p,^, and either |/i|<^^ or 0|^„.^° What these observables 
measure is correlations and their specification singles out a point in the reduced phase space. 
However, we do not have to single out a time variable to define them. The same is true 
in quantum theory. To have complete control of physics, we need to construct the physical 
Hilbert space T^phy and introduce on it a complete set of Dirac observables. Again, both 
these steps can be carried out without singling out as emergent time. For example, the 
scalar product can be constructed using group averaging which requires only the knowledge 
of the full quantum constraint and its properties, and not its decomposition into a 'time 
evolution part' 9? and an operator on the 'true degrees of freedom'. Once the scalar 
product is constructed, we can introduce a complete set of Dirac observables consisting of 
P(^, and |/i|0„ or 0|^^. Again, what matters is the correlations. This and related issues 
have been discussed exhaustively in the quantum gravity literature in relativity circles. In 
particular, a major part of a conference proceedings [6^ was devoted to it in the late eighties 
and several exhaustive reviews also appeared in the nineties (see in particular GSllGq). 



However, thanks to our knowledge of how quantum theory works in static space-times, 
singling out as the emergent time turned out to be extremely useful in practice because 
it provided guidance at several intermediate steps. In particular, it directly motivated our 
choice of L^(MBohr, -B(/i)d/iBohr) ® -^^(M, d0) as our auxihary Hilbert space; stream-lined 
the detailed definition of operators representing the Dirac observables; and facilitated the 
subsequent selection of the inner product by demanding that these be self-adjoint. More 
importantly, by enabling us to regard the constraint as an evolution equation, it transformed 
the 'frozen formalism' to a familiar language of 'evolution' and enabled us to picture and 
interpret the bounce and associated physics more easily. Indeed, following the lead of early 
LQC papers, initially we tried to use /x as time and ran in to several difficulties: specification 
of physically interesting data became non-intuitive and cumbersome; one could not imme- 
diately recognize the occurrence of the bounce; and the physics of the singularity resolution 
remained obscure. 

Our specific identification and use of emergent time differs in some respects from that 
introduced earlier in the literature. For example, in the context of the WDW theory, there 
is extensive work on isolating time in the WKB approximation (see e.g., (53l |67[). By 
contrast, a key feature of our emergent time is that it is not restricted to semi-classical 
regimes: We first isolated the scalar field as a variable which serves as a good internal 
clock away from the singularity in the classical theory, but then showed that the form of 
the quantum Hamiltonian constraint is such that can be regarded as emergent time in 
the full quantum theory, without a restriction that the states be semi-classical or stay away 



features of the LQG quantum dynamics can be traced back to the second. In particular, the bounce would 
have persisted even if we had used B(/i) in place of B{^) in the analysis presented in this paper. 
^° In the closed models, care is needed to specify the latter because they can not be defined globally. Bu t 
this issue is well- understood in the literature, especially through Rovelli's contributions 'GjI- See also 
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from the singularity. The idea of identifying an emergent time and exploiting the resulting 
'deparametrization' to select an inner product on the space of solutions to the Hamiltonian 
constraint is not new ,35^ ,5^ ,55^ ^M, ifiSi, £2i| • However, several of the concrete proposals turn 
out to have serious deficiencies (for a further discussion see, e.g., [68|,|69j). The idea of using 
a matter field to define emergent time is rather old. In the framework of geometrodynamics, 
it was carried out in detail for dust in [70|. A proposal to use a massless scalar field as time 
was also made in the framework of LQG [7l| but its implementation remained somewhat 
formal. In particular, it is unlikely that the required gauge conditions can be imposed 
globally in the phase space and the modifications in the construction of the physical scalar 
product that are necessary to accommodate a more local constructions were not spelled out. 
More recently, a massless scalar field was used as the internal clock in quantum cosmology 
in the connection dynamics framework ^] . However, the focus of discussion there is on the 
Kodama state in inflationary models. Because of the inflationary potential, the quantum 
constraint has explicit time dependence and the construction of the physical inner product 
is technically much more subtle. In particular, a viable inner product cannot depend on any 
auxiliary structure such as the choice of an 'instant of time'. These issues do not appear to 
have been fully addressed in |72j |. 



Resolutions of the Big-Bang Singularity 

The issue of obtaining singularity free cosmological models has drawn much attention 
over the years. The discovery of singularity theorems sharpened this discussion and there 
is a large body of literature on how one may violate one or more assumptions of these 
theorems, thereby escaping the big bang. Proposals include the use of matter which violate 
the standard energy conditions, addition of higher derivative corrections to the Einstein- 
Hilbert action and introduction of higher dimensional scenarios inspired by string theory. 
To facilitate comparison with the model discussed in this paper, we will restrict ourselves 
to spatially non-compact situations. 

Already in the seventies, Bekenstein investigated a model where the matter source con- 
sisted of incoherent radiation and dust, interacting with a conformal massless scalar field 
(which can have a negative energy density). He showed that Einstein's equations admit 
solutions which are free of singularities |73|. In the eighties, Narlikar and Padmanabhan 
found a singularity free solution to Einstein's equation with radiation and a negative energy 
massless scalar field (called the 'creation field') as source, and argued that the resulting 
model was consistent with the then available observations J7^. Such investigations were 
carried out entirely in the paradigm of classical relativity and the key difference from the 
standard Friedmann-Robertson- Walker models arose from the use of 'non-standard' matter 
sources. Our analysis, by contrast, uses a standard massless scalar field and every solution 
is singular in the classical theory. The singularity is resolved because of quantum effects. 

Another class of investigations starts with actions containing higher derivative terms 
which are motivated by suitable fundamental considerations. For example, to guide the 
search for an effective theory of gravity which is viable close to the Planck scale, Mukhanov 
and Brandenberger proposed an action with higher order curvature terms for which all 
isotropic cosmological solutions are non-singular, even when coupled to matter 75|. The 



modifications to Einstein's equations are thought of as representing quantum corrections. 
However one continues to work with differential equations formulated in the continuum. By 
contrast, our investigation is carried out in the framework of a genuine quantum theory 
with a physical Hilbert space, Dirac observables and detailed calculations of expectation 
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values and fluctuations. Departures from classical general relativity arise directly from the 
quantum nature of geometry. The final results are also different: While solutions in l75| 
asymptotically approach de Sitter space, in our analysis the classical big-bang is replaced 
by a quantum bounce. 

Perhaps the most well-known discussions of bounces come from the pre-big-bang cosmol- 
ogy and Ekpyrotic/Cyclic models. The pre-big-bang model uses the string dilaton action 
and exploits the scale factor duality to postulate the existence of a super-inflating pre-big- 
bang branch of the Universe, joined to the radiation dominated post-big-bang branch j^. 
However, the work was carried out in the framework of perturbative string theory and the 
transition from the pre-big-bang to post-big-branch was postulated. The initial hope was 
that non-perturbative stringy effects would enforce such a transition. However, as of now, 
such mechanisms have not been found i2l]. Although subsequent investigations have shown 
that a bounce can occur in simplified models J8[ or by using certain effective equations (see, 
e.g., [9|), it is not yet clear that this is a consequence of the fundamental theory. The 
Ekpyrotic and the more recent Cyclic models p, |^ are motivated by certain compactifica- 
tions in string theory and feature a five dimensional bulk space-time with a 4-d branes as 
boundaries. In the Ekpyrotic model, the collision between a bulk brane with a boundary 
brane is envisioned as a big bang. A key difficulty is the singularity associated with this 
collision |lO| (which can be avoided but at the cost of violating the null energy condition 
[6]). In the Cyclic model, collision occurs between the boundary branes |6|], however it has 
been shown that the singularity problems persists J76|. Thus, a common limitation of these 
models is that the branch on 'our side' of the big-bang is not joined deterministically to the 
branch on the 'other side'. In LQC by contrast, the quantum evolution is fully deterministic. 
This is possible because the approach is non-perturbative and does not require a space-time 
continuum in the background. 

Finally, the idea of a bounce has been pursued also in the context of braneworld models. 
In the original Randall- Sundrum scenario, the Friedmann equation is modified by addition 
of a term on the right side which is quadratic in density: a^/a^ = (SnG/S) p{l + p/{2a)) 
where a is brane tension. However, since a > 0, the sign of the quadratic term in p is positive 
whence a can not vanish and there is no bounce. To obtain a bounce, the correction should 
be negative, i.e., make a 'repulsive' contribution. One way to reverse the sign is to introduce 
a second time-like dimension in the bulk \JT\. However, this strategy does not appear to 
descend from fundamental considerations and the physical meaning of the second time-like 
direction is also unclear. Another avenue is to consider a bulk with a charged black hole. A 
non- vanishing charge leads to terms in the modified Friedmann equation which are negative. 
a can now vanish and a bounce can occur |73|. However it was shown that transition 
from contraction to expansion for the brane trajectory occurs in the Cauchy horizon of the 
bulk which is unstable to small excitations, thus the brane encounters singularity before 
bouncing [79|. In LQC, by contrast, while the Friedmann equation is effectively modified, 
the corrections come from quantum geometry and they are automatically negative. 

Extensions 

A major limitation of our analysis — shared by all other current investigations in quantum 
cosmology — is that the theory is not developed by a systematic truncation of full quantum 
gravity. This is inevitable because we do not have a satisfactory quantum gravity theory 
which can serve as an unambiguous starting point. The viewpoint is rather that one should 
use lessons learned from mini and midi superspace analysis to work one's way up to more 
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general situations, especially to reduce the large number of ambiguities that exist in the 
dynamics of the full theory. 

Even within quantum cosmology, our detailed analysis was restricted to a specific mini- 
superspace model. In more complicated models, differences are bound to arise. For example, 
the full solution is not likely to remain so sharply peaked on the classical trajectory till the 
bounce point and even the existence of a bounce is not a priori guaranteed, especially when 
inhomogeneities are added. However, the methods developed in the paper can be applied to 
more general situations. First, one could consider anisotropics. Now, the main structural 
difference is that the operator O will no longer be positive definite. However, a detailed 
analysis shows that what matters is just the operator |0|, obtained by projecting the action 
of O to the positive eigenspace in its spectral decomposition. Therefore, our analytical 
considerations should go through without a major modification. The numerical simulations 
will be more complicated because we have to solve a higher dimensional difference equation 
(involving 4 variables in place of 2). Another extension will involve the inclusion of non- 
trivial potentials for the scalar field. Now, generically (j) will no longer be a monotonic 
function on the classical trajectories and one would not be able to use it as 'internal time' 
globally. In the quantum theory, the operator B becomes 'time- dependent' (i.e. depends on 
0) and the mathematical analogy between the quantum constraint and the Klein-Gordon 
equation in a static space-time is no longer valid. Nonetheless, one can still use the group 
averaging procedure 35|,|5J| to construct the physical Hilbert space. For a general potential. 



a useful notion of time will naturally emerge only in the semi-classical regimes. For specific 
potentials (such as the quadratic one used in chaotic inflation) one should be able to use 
methods that have been successfully employed in the quantization of model systems J56|, |57[ 
(in particular, a pair of harmonic oscillators constrained to have a fixed total energy). 

Incorporation of spherical inhomogeneities seems to be within reach since significant 
amount of technical groundwork has already been laid J80|. Incorporation of general inho- 
mogeneities, on the other hand, will be substantially more difficult. Background dependent 
treatments have suggested that results obtained in the mini-superspace approximation may 
be qualitatively altered once field theoretical complications are unleashed (see, e.g., (8l|). 
However, already in the anisotropic case, there is a qualitative difference between pertur- 
bative and non-perturbative treatments. Specifically, if anisotropics are treated as pertur- 
bations of a background isotropic model, the big-ban g si ngularity is not resolved while if 
one treats the whole problem non-perturbatively, it is J82]. Therefore definitive conclusions 
can not be reached until detailed calculations have been performed in inhomogeneous mod- 
els. However, if a quantum bounce does generically replace the big bang singularity, it 
would be possible to explore the relation between the effective descriptions of LQG and the 
Hartle-Hawking 'no boundary' proposal J83|. For, in the effective description, the extrinsic 
curvature would vanish at the bounce. Therefore generically it may be possible to attach 
to the Lorentzian, post-bounce effective solution representing the universe at late times, a 
Riemannian pre-bounce solution without boundary. If so, it would be very interesting to 
analyze the sense in which this Riemannian solution captures the physics of the pre-bounce 
branch of the full quantum evolution. 

Finally, it is instructive to recall the situation with singularities in classical general rela- 
tivity. There, singularities first appeared in highly symmetric situations. For a number of 
years, arguments were advanced that this is an artifact of symmetry reduction and generic, 
non-symmetric solutions will be qualitatively different. However, singularity theorems by 
Penrose, Hawking, Geroch and others showed that this is not correct. An astute use of the 
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differential geometric Raychaudhuri equation revealed that singularities first discovered in 
the simple, symmetric solutions are in fact a generic feature of classical general relativity. 
A fascinating question is whether the singularity resolution due to quantum geometry is 
also generic in an appropriate sense J8J]. Is there a general equation in quantum geometry 
which implies that gravity effectively becomes repulsive near generic space-like singulari- 
ties, thereby halting the classical collapse? If so, one could construct robust arguments, 
now establishing general 'singularity resolution theorems' for broad classes of situations in 
quantum gravity, without having to analyze models, one at a time. 

Acknowledgments: We would like to thank Martin Bojowald, Jim Hartle and Pablo 
Laguna for discussions. This work was supported in part by the NSF grants PHY-0354932 
and PIIY-0456913, the Alexander von Humboldt Foundation, and the Eberly research funds 
of Penn State. 



APPENDIX A: HEURISTICS 

Quantum corrections to the classical equations can be calculated using ideas from a 
geometric formulation of quantum mechanics where the Hilbert space is regarded as (an 
infinite dimensional) phase space, the syrnplectic structure being given by the imaginary part 
of the Hermitian inner product (see, e.g., |85|). This 'quantum phase space' has the structure 
of a bundle with the classical phase space as the base space, and all states with the same 
expectation values for the canonically conjugate operators (g*, pi) as an (infinite dimensional) 
fiber. Thus, any horizontal section provides an embedding of the classical phase space into 
the 'quantum phase space'. In the case of a harmonic oscillator (or free quantum fields) 
coherent states constitute horizontal sections which are furthermore preserved by the full 
quantum dynamics. In the semi-classical sector defined by these coherent states, the effective 
Hamiltonian coincides with the classical Hamiltonian and there are no quantum corrections 
to classical dynamics. For more general systems, using suitable semi-classical states one 
may be able to find horizontal sections which are preserved by the quantum Hamiltonian 
fiow to a desired accuracy (e.g. in an h expansion). The effective Hamiltonian governing 
this flow — the expectation value of the quantum Hamiltonian operator in the chosen states, 
calculated to the desired accuracy — is generally different from the classical Hamiltonian. 
In this case, dynamics generated by the effective Hamiltonian provides systematic quantum 
corrections to the classical dynamics 25| (see also |26l. I27ll34l|). 



This procedure has been explicitly carried out in LQC for various matter sources j2 
86|. For a massless scalar fleld, the leadingorder quantum corrections are captured in the 



following effective Hamiltonian constraint |86l |: 

Ceff = ^bl^ SlYi\^loC)+%'KGB{p)pl (Al) 

where B{p) is the eigenvalue of l/|pp/^ operator given by ()4.3|) . ^^ For |yu| ^ /i^, B{p) can 



^^ In the literature, eigenvalues B[p) often contain a half-integer j, a parameter representing a quantization 
ambiguity. In view of the general consistency arguments advanced in [Sg, we have set its value to its 
minimum, i.e. j = 1/2. 
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be approximated as 
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The leading order term is l/p^^"^, thus B(p) quickly approaches its classical value for \fi\ ^ /io, 
corrections being significant only in the 'genuinely quantum region' in the vicinity of /x = 0. 

From now on, we will ignore the quantum corrections to B{p). 
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FIG. 17: Expectation values (and dispersions) of |/i|(^ are plotted near the bounce point, together 
with with classical and effective trajectories (fainter and darker dots, respectively). While the 
classical description fails in this region, the effective description provides an excellent approximation 
to the exact quantum evolution. In this plot, p^ = 3000 and e = 2fio- 



To obtain the equations of motion, we need the effective Hamiltonian Ties- As usual it is 
obtained simply by a rescaling of Ccs which gives TCes the dimensions of energy and ensures 
that the matter contribution to it is the standard matter Hamiltonian: 
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Then, the Hamilton's equation for p become: 

Sn^GdHes 2\p\y' 
p = [p, Hcs} = 5 7^— = 

6 OC 7/io 

Further, since the Hamiltonian constraint implies that Tics of ()A3|) vanishes, we have: 
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which, on using Eq. ()A4|) . provides the modified Friedmann equation for the Hubble param- 
eter H: 

H =Tl = ^rP{^- 1^ where Pent = 0^22 — • ^^ 

To obtain the dynamical trajectory, we also need the Hamilton's equation for (p, 

0={</',K4 = ^- (AT) 

By combining Eq. ()A4|) and ()A7|) we obtain the effective equation of motion in the /i-0 
plane: 

The classical Friedmann dynamics results if we set pcrit = 00. Eqs. ()A6j) and ()A8jl suggest 
that the LQC effects significantly modify the Friedmann dynamics once the matter density 
reaches a critical value, Pcrit- In the classical dynamics, the Hubble parameter H can not 
vanish (except in the trivial case with p^ = 0). In the modified dynamics, on the other hand, 
H vanishes at p = pcrit- At this point, the Universe bounces. Thus, the bounce predicted 
by Eq. ()A6|) has its origin in quantum geometry. (The critical value /x at which this bounce 
occurs is given by /icrit = {^/Q/^'nGffl) p^ /ig.) As pointed out in the main text, a physical 
limitation of the present framework is that if p^ is chosen to be sufficiently large, the critical 
density pcrit can be small. 

To obtain the effective equations, several approximations were made |25l . ISq which are 
violated in the deep Planck regime. Nonetheless, the resulting picture of the bounce is con- 
sistent with the detailed numerical analysis. In fact, within numerical errors the trajectory 



P(0) = ^ I exp I \/^(0 - 0o) I + D{p^) exp I -7^(0 - 0o) 1 I (A9) 




obtained by integrating Eq ()A8|) approximates the expectation values of \fi\^ quite well. (As 
in the main text, D{p^) = p^^pj,/12TiGfi^). An illustrative plot of this generic behavior is 
shown in Fig. El Therefore, in retrospect, this analysis can be taken as a justification for 
the validity of the approximation throughout the evolutionary history of semi-classical states 
used in this paper. However, by its very nature, the effective description can not reproduce 
the interesting features exhibited by quantum states captured in Figs. IT^HTHl 

However, the effective description can be used to provide an intuitive understanding of 
the behavior of various uncertainties discovered through numerical analysis. Let us first note 
that the position of the bounce point depends linearly on the value oip^. Next, consider two 
near-by solutions with slightly different p^ which asymptote to the same WDW solution for 
the expanding branch in the distant future. We wish to know the way in which {6 fi / fi) {(p) := 
((/ii — /i2)//i2)(0) changes in the backward evolution as the two wave functions asymptote 
to WDW solutions in the distant past. This relative difference can be found using Eq ()A9|) 
and is given by 

^ = 2^+f^V, (AlO) 

P P^ \P<P J 
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FIG. 18: The darker 'lattice' show the difference between final and initial relative spreads A///// 
for states obtained by evolving initial data of Method II. The upper, fainter 'lattice' shows the 
heuristic bound 6n/iJ, given by Eq ()A10|) . 



where 5p^ is the difference between values oi p^ of the two classical trajectories. A heuristic 
estimate on the relative difference in /x can be compared with the relative dispersion Afi/n 
obtained from the Method II in the (p evolution of section IV Bl It turns out that estimate 
in Eq ()A10|) provides a reasonably good upper bound to the relative dispersion found nu- 
merically (see Fig. ITH|) . A similar comparison can be made for Method III of section IV Bl 
In this case the corresponding solution to Eq ()A8j) is 
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A comparison with the relative dispersions in numerical analysis is shown in Fig. ^1 As in 
the case of construction of coherent states via Method II, above estimate serves as an upper 
bound for relative dispersions computed by numerical analysis. 

Finally, in the numerical analysis an important issue concerns with the behavior of disper- 
sions of our Dirac observables /i^ and p<^ and the product AcpAp^. Intuitive understanding 
of our numerical results of Fig. ^J can be gained by casting Eq. ()A8|) in the form 
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FIG. 19: The darker 'lattice' show the difference between final and initial relative spreads A///// 
for states obtained by evolving initial data of Method III. The upper, fainter 'lattice' shows the 
heuristic bound Suf/j, given by Eq (|A12() . 



Since Afi/fi can be determined numerically, we can then estimate A(j) throughout the evolu- 

/ ^ \-l/2 

tion. The factor 1 — is approximately equal to unity for p <^ pcrit or equivalently 

y Pcrit / 

/ \-l/2 

for u ^ yUcrit- Howcvcr near the bounce point, 1 — ^ 1. In Method II and 

\ pcrit / 

Method III of constructing initial states described in Sec. IV B 21 this change compensates 
the corresponding decrease in A/i//i and leads to a nearly constant value of A0. However, 
for Method I, since A/i//i increases monotonically, the fluctuation A0 increases significantly 
near the bounce point. 



APPENDIX B: ISSUES IN NUMERICAL ANALYSIS 



In this Appendix we will spell out the way in which the WDW limit of the eigenfunctions 
of G were found. 

Consider a general eigenfunction e^^ of ( for cu'^ > vrG/S). It is always a linear com- 
bination of basis functions e|;i.|,e_|;i.| (where k"^ = 3/(167rG)ti;^ — 1/16) defined in Eq. (j3.5|) . 
For later convenience, let us express the linear combination as: 



e^ = r+e'(^+")ei,i + r-e'(^-")e_ 



■£|fc| -r ■' c ■ -c.ifci , (Bl) 

where r±,a, P are real numbers. Since each e_|_|j;.| is a product of an 'amplitude' \fi\^ /An and 
a 'phase' e=''*l'^l''^l'^l, it is natural to rescale e^: 



L(/^) :=47r|/i| 46^ 
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In terms of coefficients defined in ()B1|) . we have 

e^(/i) = e*^ ((r"'' + r~) cos(q; + |/i;| In |/i|) + i(r'^ — r^) sin(a + |A;| In |/i|)) . (B3) 

Tfie values of e^(/x) trace out an ellipse on the complex plane, parameterized by In |/i|. The 
length of semi-major and semi- minor axis of this ellipse is equal to, respectively 

r^ -\- r~ = sup |e^| \r^ — r^\ = inf |e^^| , 



M 



A* 



whereas the phase a is related to positions of maxima of e^ as follows 

|e^| = r+-|-r^ -v^ a + \k\ln\fi\ = rnr , n G Z . (B5) 

The remaining phase f3 is just the phase of e^ at maximum. The sign of (r"*" — r~) is, on the 
other hand, determined by direction of the rotation of the curve as /i increases. 

The method specified above allows us to calculate the decomposition in e^ basis of a 
function e^^ specified in the form of numerical data (i.e. array of values at sufficiently large 
domain). The same algorithm can be applied to identify the WDW limit of any eigenfunction 
of the LQC operator 6. Indeed given an eigenfunction e^{fi) supported on the lattice C\^\ 
(or >C-|e|) one can again define e^^ analogously to e^^ and find its (local) extrema for large 
/i. (For definiteness, we restrict our consideration to finding the limit on the positive n 
side, however this method can be used of course also for the negative /i domain.) Next, 
the positions of extremas and values of e,^ at them can be used to calculate coefficients 
r"*" ± r~,a,/3 at each extremum independently. If e^^ — >■ e^ then these coefficients form 
sequences ({r+ + r~}j, {r+ — r~}j, {a}i, {P}i) which converge to the analogous coefficients 
corresponding to e^ a.s fi —>■ oo. Finding the WDW limit of e^j reduces then to finding the 
limit of ({r+ + r-}j, {r+ - r"}i, {a}i, {P}i). 

In actual numerical work the following method was used: 

• After given eigenfunction e^^ was calculated using ()4.(jj] the positions of extrema {/i}j 
were found. 

• Around the extrema the function e^ supported on C\s\ (or C-\e\) was extended to 
neighborhoods of {/x}, via polynomial interpolation. Then the positions and values 
of extrema were recalculated with use of this extension. This allowed us to construct 
sequences converging to the WDW limit much more quickly than the ones constructed 
in the first step. The motivation for this construction is the expectation that for 
sufficiently large /i the values of e^^ at C\e\ should be good estimates of its WDW limit 
(being regular function defined on entire IR+). 

• Extrema found in previous step were next used to calculate sequences ({r++r^}j, {r+— 
r~}i, {a}i, {/9}j) in a way analogous to that specified by Eqs. ()B4|) and ()B5|) and in the 
description below them. 

• Finally the limits of coefficients at l//i — ^ were calculated by polynomial extrapola- 
tion. 
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APPENDIX C: AN ALTERNATE PHYSICAL HILBERT SPACE 



In this Appendix we will construct a physical Hilbert space H'^^^n LQC which is qual- 
itatively different from the spaces Hp^y constructed in section II V Bl In its features, it 
interpolates between these and the Hilbert space Tip^j" of section IIIIBI For completeness, 
we will first explain why a new representation of the algebra of Dirac observables can arise, 
then summarize the results and finally compare and contrast them with those obtained 
in sections IIIIBI and IIV Bl The first part is somewhat technical but we have organized the 
presentation such that the readers can go directly to the summary without loss of continuity. 

Let us begin by recalling the situation for a general system with a single constraint C. In 
the refined version of Dirac quantization J36|, one introduces an auxiliary Hilbert space T^aux, 
and represents the constraint by a self-adjoint operator C on it. The technically difficult 
task is to chooses a dense sub-space $ of Tiaux such that for all /,(?€$, 



(^J 



dX{e-'''-f\ 



(CI) 



is a well-defined element of $*, such that the action {"^flg) of (\E'/| G $* on \g) G $ yields a 
Hermitian scalar product on the space of solutions (\I'/| to the quantum constraint (see, e.g. 
35|, ISy, |87| ) . Results in section ITV Bl were obtained using L^(RBohr, B{^) dyUBohr) ®i^^(M, d0) 
for Haux, and the space of rapidly decreasing functions /(/i, (p) in this Tiaux for $. In the 
LQC literature, $ is sometimes called Cyl and $* is taken to be its algebraic dual, denoted 
by Cyr. 

The construction given above is rather general. For the model under consideration, we 
can extend this construction by using an entirely different subspace of $* for the auxiliary 
Hilbert space. This is possible because by duality the action of C can be extended to all of 
$*. Let us set T^aux •= L'^i^'^, -B(/x)d/id0). This is a subspace of $* because each \E' G Ti.'^^ 
defines a linear map from $ to C: 



/oo 
d05(/x)^(M,0)/(/x,0) Vv&gK, 



(C2) 



where the sum over fx converges because / G $ has support only on a countable number of 
points on the /i-axis and a rapid fall-off. The dual action of C on T^aux can now be calculated: 
Since 
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(C3) 



it follows from the definitions of C^(/i) that 
52^ 



C^){fx,(t>) = ^ + [Bifi)]-' (C+(^)vl>(/i + 4/i„0) + C°(/i)vl>(;x,0) + C-(/i)vl>(^-4/i„)) 
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It is straightforward to verify that G and C are self-adjoint on Ti^ux- Therefore, we can 
carry out group averaging on Ti^ux ^'^d obtain a new physical Hilbert space. As in the 
WDW theory of section IIIIB[ there are two superselected sectors. We will work with the 
positive frequency sector and denote it by Hp^y The Dirac observables p^ and l/xj^^ on $ 
act by duality on Ti^ux ^'^d descend naturally to 'H'^-^y 

The final results can be summarized as follows. The new physical Hilbert space '^^phy i^ 
the space of functions \l'(yU,0) satisfying the positive frequency equation: —id^'^ = -\/9^, 
with finite norm: 

\\n%y= [ df,B{f,)\^{f,,<P)\' (C5) 

J 4>=<t>o 
and the action of the Dirac observables is the standard one: 

i;:j;„v^(/x,0) = e^-^(<^-*°)Hvl/(;.,0„), and p^^i;(^,c^) = -^h^^^^. (C6) 

Note that the final physical theory is different from both the WDW theory of 1111 Bl and 
the 'standard' LQC theory of llVBl Since the inner product ()C5|) involves an integral rather 
than a sum, states \l/(/U, 0) now have support on continuous intervals of the /x-axis as in the 
WDW theory, rather than on a countable number of points as in LQC. However, the states 
satisfy the LQC type positive frequency equation —id^'^ = \/9^, where the operator on 
the right is the square-root of a positive, self-adjoint difference operator rather than of a 
differential operator, and the measure determining the inner product also involves -B(/i) from 
LQC rather than B(yu) from the WDW theory. Thus, the dynamical operator is the same as 
in LQC. In particular, as in section\^ the quantum states exhibit a big bounce. However, 

since typical states are continuous on the /i-axis, the spectrum of the Dirac observable \fi\ct,^ 
is now continuous. In essence, the states \E'(/i, 0) have support on 2-d continuous regions 
of the /i — plane as in the WDW theory but their dynamics is dictated by a difference 
operator as in LQC. When the cosmological constant is non-zero, the analog of this physical 
Hilbert space appears to be a natural home to analyze the role of the Kodama state in 
quantum cosmology |72l |. 

In the literature on 'polymer representations', non-relativistic quantum mechanics of 
point particles and the quantum theory of a Maxwell theory have been discussed in some 
detail|48|, |89|. In the first case, the standard Schrodinger Hilbert space L^(R, dx), and in 
the second case, the standard Fock space turned out to be subspaces of Cyl* which were 
especially helpful for semi-classical analysis. The present auxiliary Hilbert space Tiaux '^ ^* 
is completely analogous to these. Therefore the resulting Hp^y may be more useful for 
semi-classical considerations. Indeed, since it does not refer to any e, no coarse graining is 
required to carry out the semi-classical analysis. Therefore the analog of this construction 
may well be useful in more general contexts in full LQG. 
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